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ABSTRACT 

It is now well established that all galaxies with a massive bulge component harbour a cen- 
tral supermassive black hole (SMBH). The mass of the SMBH correlates with bulge prop- 
erties such as the bulge mass and the velocity dispersion, which implies that the bulge and 
the central black hole of a galaxy have grown together during the formation process. As 
part of an investigation of the dependence of the SMBH mass on bulge types and forma- 
tion mechanisms, we present measurements of SMBH masses in two pseudobulge galax- 
ies. The spiral galaxy NGC 3368 is double-barred and hosts a large pseudobulge with a tiny 
classical bulge component at the very centre. The SO galaxy NGC 3489 has only a weak 
large-scale bar, a small pseudobulge and a small classical bulge. Both galaxies show weak 
nuclear activity in the optical, indicative of the presence of a supermassive black hole. We 
present high resolution, adaptive-optics-assisted, near-infrared integral field data of these two 
galaxies, taken with SINFONI at the Very Large Telescope, and use axisymmetric orbit mod- 
els to determine the masses of the SMBHs. The SMBH mass of NGC 3368, averaged over 
the four quadrants, is <M.> - 7.5 x 10^ Mq with an error of 1.5 x 10^ Mq, which mostly 
comes from the non-axisymmetry in the data. For NGC 3489, a solution without black hole 
cannot be excluded when modelling the SINFONI data alone, but can be clearly ruled out 
when modelling a combination of SINFONI, OASIS and SAURON data, for which we obtain 
M. = (6.00;^Q 54lstat ± 0.64|sy.s) X 10'' Mq. Although both galaxies seem to be consistent with 
the M,-<j relation, at face value they do not agree with the relation between bulge magni- 
tude and black hole mass when the total bulge magnitude (i.e., including both classical bulge 
and pseudobulge) is considered; the agreement is better when only the small classical bulge 
components are considered. However, taking into account the ageing of the stellar population 
could change this conclusion. 

Key words: galaxies: kinematics and dynamics — galaxies: bulges — galaxies: individual 
(NGC 3368, NGC 3489) 



1 INTRODUCTION 

Bulges located in the central regions of disc galaxies are commonly 
identified as the region where the excess light above the outer expo- 
nential disc dominates the surface brightness profile. Bulges were 
generally regarded as scaled-down versions of elliptical galaxies, 
probably formed via minor galaxy mergers. There is now evidence 
that there is also a second type of central structure, the so-called 

* Based on observations at the European Southern Observatory VLT 
(078.B-0103(A)), and on service observations made with the William Her- 
schel Telescope operated on the island of La Palma by the Isaac Newton 
Group in the Observatorio del Roque de los Muchachos of the histituto de 
Astrofi'sica de Canarias. 
t E-mail: nnowak@mpe.mpg.de 



pseudobulges, which are more similar to mini-discs than to mini- 
ellipticals. They were first introduced by ,Kormendy| ^1993^) , and 
[Kormendy & Kennicutt, Jn| j2004[ > review properties and forma- 
tion mechanisms and present a number of examples. Pseudobulges 
are thought to be the result of secular evolution and can be identi- 
fied e.g. through the presence of disc-like structure (nuclear spirals, 
bars or rings), flattening similar to that of the outer disc, rotation- 
dominated kinematics, exponential surface brightness profiles or 
young stellar populations. As the formation mechanisms of classi- 
cal and pseudobulges are fundamentally different and can happen 
independently, galaxies could harbour both types of bulges ( Erwm| 
let aL|2003| |Athanassoula|2005[ [Erwin|2008^ . But this fundamen- 
tal difference between the formation mechanisms also leads to the 
question whether and how a central black hole grows inside a pseu- 
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Table 1. Properties of the galaxies NGC 3368 and NGC 3489. As both galaxies host a composite bulge, the effective radius 
and the K-banA bulge magnitude Mic are given for the photometric bulge and the classical bulge component. 
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dobulge and how the mass of the black hole relates to pseudobulge 
properties. Supermassive black holes (SMBHs) in elliptical galax- 
ies and classical bulges are known to follow tight correlations with 
luminosity (e.g. Kormendy & Richstone||1995] [Marconi & Hunt| 
[2003l >, mass (Haring & Rix 2004) and velocity dispersion (M.-cr 
relation, [Gebhardt et al. 2000b, Ferrarese & Merritt„2000 l of the 
bulge. It is not clear whether pseudobulges follow the same rela- 
tions, for several reasons: 1. There are only very few direct SMBH 
mass measurements in pseudobulges. 2. The concept of pseudob- 
ulges is relatively new and the classification criteria therefore dif- 
fer somewhat from author to author. 3. The fact that at least some 
galaxies could contain both bulge types (composite bulges) makes 
the classification and correlation studies even more complicated. 
The composite bulges need to be decomposed properly in order to 
find out with which property of which bulge component the SMBH 
mass correlates. [Kormendy I ( [200 1 ) did not find any dependence of 
the M.-cr relation on the mechanism that feeds the black hole. In 
contrast ' Hu[ ( |2008| l finds that the black holes in pseudobulges have 
systematically lower masses than black holes in classical bulges 
and ellipticals with the same velocity dispersion. Both studies suf- 
fer from small number statistics and unclear classification issues. 
For low-mass galaxies without classical bulge (i.e. likely hosts of 
a pseudobulge) and with virial SMBH mass estimates [Greene, Ho[ 
[& Barth[ ( 12008) found no deviation from the M.-cr relation, but a 
likely disagreement with the M.-Mbuigc relation. [Gadotti & Kauff-[ 
[mann[ ( [2b09^ found for a large number of SDSS galaxies that pseu- 
dobulges, classical bulges and ellipticals cannot follow both the 
M.-a and the M.-Mbuige relation at the same time. As they esti- 
mated M. from the M.-Mbuige relation by [Haring & Rix[ l [2004^ it 
is not clear whether their pseudobulges follow a different M.-cr or 
M.-Mbiiigc relation or both. 

In this paper we present a thorough analysis and derivation 
of the black-hole masses via extensive stellar dynamical modelling 
of NGC 3368, a double-barred spiral galaxy ( Erwin|2004^ of type 
SAB(rs)ab with a well-defined pseudobulge and a very small clas- 
sical bulge component, and NGC 3489, a SAB(rs)0-l- galaxy with 
a weak large bar, a small pseudobulge and a similar-sized classi- 
cal bulge. All important parameters of the two galaxies are listed 
in Table ^ Using high-resolution imaging we are able to iden- 
tify and decompose the pseudobulge and classical bulge compo- 
nents. High-resolution adaptive-optics assisted near-IR integral- 
field spectroscopy enables us to model each quadrant separately. In 
contrast to our two previous studies of elliptical galaxies ( |Nowak[ 
[et al. 2007 2008), non-axisymmetries may play a larger role due to 
the barred nature of the galaxies. 

The nucleus of NGC 3368 is weakly active and can be classi- 
fied as a LINER2 based on optical emission line ratios ( Ho, Filip^ 
[penko & Sargent|I997^ . [Maoz et al.[ l |2005) and [Maoz[ ( [2()07[ ) report 
long-term UV variations, which suggests the presence of an AGN 



and thus a SMBH. NGC 3489 has a weak LINER/HII transition 
type or Seyfert 2 nucleus ( |Ho et al.|1997| ). 

We adopt a distance to NGC 3368 of 10.4 Mpc throughout the 
paper based on surface brightness fluctuation measurements ( |Tonry[ 
[et al.|20 0r). At this distance, I arcsec corresponds to ~ 50 pc. For 
NGC 3489 we adopt a distance of 12.1 Mpc (~ 59 pc arcsec"'), 
also based on the measurements of jTonry et al.[pOOT) . If the M.- 
cr relation jTremaine et al.|2002| > applies, both galaxies are close 
enough to resolve the sphere of influence of the black hole from 
the ground with adaptive optics. 

This paper is organised as follows: In §2 we discuss the mor- 
phology of the two galaxies, including photometric (and spectro- 
scopic) evidence for pseudobulges and classical bulges. The spec- 
troscopic data, including stellar kinematics, gas kinematics and line 
strength indices are described in §3. The stellar dynamical mod- 
elling procedure and the results for the SMBH mass of each galaxy 
are presented in §4 and §5. §6 summarises and discusses the results. 



2 IMAGING 

2.1 Identifying Discs, Pseudobulges, and Classical Bulges 

As mentioned in the Introduction, both NGC 3368 and NGC 3489 
have complex morphologies; they are not simply an exponential 
disc plus a classical bulge. In this section, we explain our approach 
for analysing the morphology of these galaxies, and how we de- 
compose their surface brightness profiles for purposes of dynamical 
modelling. 

Since discs and bulges can have different mass-to-light ratios 
T, we need to separate these components for modelling purposes. 
But we also need to ensure that what we call the "bulge" really is 
distinct from the disc, and not simply a higher-surface-brightness 
extension of the disc. NGC 3368 has already been classified by 
[Drory & Fi'sherl pOOTl l as a galaxy where the bulge is a "pseu- 
dobulge", so we need to consider whether this galaxy even has a 
distinct bulge; if it does have a pseudobulge, we need to identify 
that. 

Our overall approach is as follows: First, we perform a "naive" 
global bulge/disc decomposition, where we treat the entire surface- 
brightness profile as the combination of an outer exponential (the 
"main disc") and a Sersic profile. The latter component is the "pho- 
tometric bulge," which would traditionally be considered the bulge 
- i.e., the kinematically hot spheroid - of the galaxy. Second, we 
focus on the photometric bulge region (that is, where the Sersic 
component dominates the modelled global light profile), and ex- 
amine the morphology and stellar kinematics. In both galaxies, we 
find evidence that these regions are predominantly disclike and thus 
not classical bulges. This argues that we should treat most of the 
photometric bulge region as part of the disc. Finally, we also find 
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Figure 1. Isophotal maps of NGC 3368 on different scales, all with logarithmic intensity scaling, (a) SDSS r-band isophotes, showing the outer disc. The gray 
box outlines the region shown in the next panel, (b) A'-band isophotes from the image of |Knapen et al.|(2003 ^, showing the outer bar and the lens just outside, 
(c) A'-band isophotes, showing the region inside the outer bar, including the bright inner pseudobulge region, (d) A'-band isophotes, showing the elliptical 
contours of the inner pseudobulge, the inner bar, and the central bulge inside. The dashed blue square shows the approximate field of view of our lOOmas 
SINFONI observations. For reference, we indicate the adopted major-axis position angle (172°) with the diagonal red Une in panels b-d. 



evidence that the central few hundred parsecs contain an additional 
component: a central light excess above the disclike bulge region 
associated with kinematically hot stellar kinematics and rounder 
isophotes. It is this last component which we call the "classical 
bulge," and which we treat as a separate stellar component in our 
modelling. 

We thus consider both galaxies to be similar to NGC 2787 
and NGC 3945, where ,Erwin et al.| pOOB] ! showed that the pho- 
tometrically defined "bulges" of both galaxies were really bright 
"inner discs," with much smaller (and rounder) bulges inside the 
inner discs. Erwin (2008 ) and Erwin et al. (2009, in prep) discuss 
more examples of such "composite bulge" galaxies, andjAthanas-| 



|soula| ( [2005( 1 provides a theoretical context for such systems. 

To summarise: Both galaxies appear to consist of a disc with 
an outer exponential profil^and a steeper inner profile, along with 



' This is a simplification, since the profile of the disc at very large radii 
may change jErwin et al.|2008| , but this has no eft'ect on our modelling. 



a small central excess with rounder isophotes and stellar kinematics 
which we consider to be the (classical) bulge. The steep inner part 
of the disc can be considered a discy pseudobulge, but for mod- 
elling purposes we treat it as just the inner part of the disc. 

2.2 Imaging Data and Calibrations 

The imaging data we use comes from a variety of sources, including 
the Two-Micron All-Sky Survey (2M ASS;|Skrutskie et al.|2006l ), 
the Sloan Digital Sky Survey (SDSS; !^^: et al.||2000[ l, and the 
Hubble Space Telescope archive. We also use near-IR images taken 
with the Isaac Newton Group Red Imaging Device (INGRID, a 
1024^ near-IR imager with 0.24 arcsec pixels) on the William Her- 
schel Telescope: a A^-band image of NGC 3368 from |Knapen et al.| 
(2003), available via NED, and an //-band image of NGC 3489 ob- 
tained during service/queue time (February 11, 2003). The seeing 
for the INGRID images was 0.77 arcsec FWHM for NGC 3368 and 
0.74 arcsec FWHM for NGC 3489. Finally, we also use AT-band im- 
ages created from our SINFONI datacubes (see Fig.[T3](. 
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Figure 2. Isophotal ellipse fits for NGC 3368, plotted on logarithmic (left panel) and linear (right panel) scales. Black points are for the SDSS r-band image, 
red points are for the X^-band image of Knapen et al. ( 2003 1, purple points are for the NICMOS2 F160W image and green points are for the SINFONI image 
(derived from our datacube). Major morphological features are marked ("Bar 1" = outer bar, "Bar 2" = inner bar). The dashed line in the top panels indicates 
our adopted position angle (172°) for the galaxy disc. 



The SDSS r-band images are used for measuring the shape of 
the outer-disc isophotes, which helps us determine the most likely 
incUnation for each galaxy. The innermost isophote shapes and sur- 
face brightness profiles are determined from the HST archival im- 
ages. For NGC 3368, we choose a NICMOS2 F160W image (PI 
Mulchaey, proposal ID 7330) for that purpose; with F450W and 
F814W WFPC2 images (PI Smartt, proposal ID 9042) we con- 
struct high-resolution colour maps and attempt to correct the NIC- 
M0S2 image for dust extinction. For NGC 3489, we used F555W 
and F814W WFPC2 images (PI Phillips, proposal ID 5999). We 
attempted to correct the NICMOS2 F160W image (for NGC 3368) 
and the WFPC2 F814W image (for NGC 3489) for dust extinction, 
following the approach of [Barth et al.| ( |2001^ and [Carollo et al.| 
( |I997^ . This involved creating a V - H colourmap for NGC 3368 
and a V - / colourmap for NGC 3489, then generating correspond- 
ing Am and A/ extinction maps and correcting the NICM0S2 and 
WFPC2 F814W images. The results were reasonably successful 
for NGC 3489, but less so for NGC 3368, perhaps due to the much 
stronger extinction in the latter galaxy. 



The 2MASS images are used primarily to calibrate the IN- 
GRID near-IR images. Since the latter suffer from residual sky- 
subtraction problems, we calibrate them by matching surface- 
brightness profiles from the INGRID images with profiles from the 
appropriate 2MASS images (A'-band for NGC 3368, //-band for 
NGC 3489), varying both the scaling and a constant background 
ofi'set until the diff'erences between the two profiles are minimized. 
We then carry over this calibration to surface-brightness profiles 
from the HST images: i.e., we calibrate the NICMOS2 F160W 
profile to AT-band for NGC 3368 by matching it to the (calibrated) 
INGRID AT-band profile (including a background offset), and simi- 
larly match the WFPC2 F814W profile to INGRID //-band profile 
for NGC 3489. Profiles from the SINFONI A:-band images are then 
calibrated by matching them to the appropriate calibrated HST pro- 
files. 



2.3 NGC 3368 

2.3.1 Morphological overview 

NGC 3368 is a relatively complex spiral galaxy, with a number of 
different stellar components. |Erwin| ( |2004[ l argued that the central 
regions of NGC 3368 included at least three distinct components: 
an outer bar with semi-major axis a a 61-75 arcsec (4.4-5.4 kpc, 
deprojected), an "inner disc" extending to a k 21-30 arcsec (1.1- 
1.6 kpc, deprojected), and an inner bar with a x 3.4-5.0 arcsec 
(200-300 pc, deprojected). As noted above, this set of nested struc- 
tures is very similar to that of the double-barred galaxy NGC 3945 
l |Erwin & Sparke|1999^ , where [Erwih et al.|p003l > found that the 
galaxy's "photometric" bulge could be decomposed into a bright, 
kinematically cool disc (first noted by Kormendy 1982 1 with an ex- 
ponential profile and a much smaller, rounder object dominating 
the inner few hundred parsecs - apparently a central, spheroidal 
bulge. 

The isophotes of NGC 3368 are shown in Fig.[T]for diflrerent 
scales. The isophotal ellipse fits to ground-based and HST near-IR 
images are shown in Fig.|2] In the inner region (r ~ 3-4 arcsec), the 
ellipticity rises to a local maximum of ~ 0.3 and the isophotes are 
closely aligned with the outer disc. Even further in, inside the in- 
ner bar (semi-major axis a <2 arcsec), the isophotes become quite 
round, with a mean ellipticity of x 0.1. The isophotes in this re- 
gion also twist significantly; inspection of both the NICM0S2 im- 
age and our SINFONI datacubes indicate that this twisting is pro- 
duced by strong dust lanes on either side of the galaxy centre. The 
true (unextincted) ellipticity in this region is probably close to 0. 
This suggest that NGC 3368 harbours a small classical bulge, again 
in analogy to NGC 3945. The size of the classical bulge region is 
much larger than the NICMOS PSF and is therefore well resolved 
and an AGN can thus be excluded. In addition we do not find emis- 
sion lines characteristic for an AGN in our SINFONI spectra (see 
below). 
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Figure 3. Global bulge-disc decomposition of NGC 3368. The data points 
(circles) are the major-axis /T-band profile, combining the HST-N\CM0S2 
F160W image (r < 7.1 arcsec) with the ^T-band image of [Knapen et al.| 
j2003) for 7.2 < r < 86 arcsec and the SDSS r-band image for r > 86 arc- 
sec, all calibrated to K. Also shown is the best Sersic + exponential fit to 
the data and the residuals (bottom panel). The Sersic component represents 
the "photometric bulge," which dominates the light at r < 50 arcsec. (Most 
of this, we argue, is really a part of the disc; see text for details.) 
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Figure 5. Local estimates of the ratio of ordered motion to random mo- 
tion in stellar kinematics along the major axis of NGC 3368, within the 
photometric bulge region {r < 50 arcsec), based on the long-slit data of 
[Heraudeau et al.|jl999^ . We first deproject the stellar velocities to their in- 
plane values (correcting from the observed PA of 5° to our adopted major- 
axis PA of 172°), then divide them by the observed velocity dispersion val- 
ues. Ratios < 1 indicate a dominance of velocity dispersion over bulk ro- 
tation; ratios > 1 indicate kinematically cooler regions. The vertical green 
dotted lines indicate the seeing of the original spectroscopic observations 
(FWHM= 1.8 arcsec). 
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Figure 4. Bulge-disc decomposition of the inner photometric bulge of 
NGC 3368. The data points (circles) are the major-axis A'-band profile for 
the inner r < 8.5 arcsec, combining our SINFONI data cube (r < 0.3 arc- 
sec) with the HST-mCMOS2 F160W image, both calibrated to K. Also 
shown is the best Sersic + exponential fit to the data at r < 8 arcsec and 
the residuals (bottom panel); in this fit, the Sersic component represents the 
"classical bulge,", while the exponential is the inner part of the photometric 
bulge - that is, it is the steep inner part of the galaxy disc. 



2.3.2 Bulge-disc decomposition 

Fig. |3] shows a global bulge-disc decomposition for NGC 3368, in 
which the photometric bulge (the Sersic component of the fit) dom- 
inates the hght at r < 50 arcsec. 

Fig. |4] shows a Sersic -I- exponential decomposition of the in- 
ner r < 10 arcsec region. In this decomposition, we are now treat- 
ing what we previously identified as the photometric bulge as a 



disc-like component (i.e. the pseudobulge), which has (compared 
to the outer disc) a relatively steep exponential profile, plus a cen- 
tral Sersic excess. Note that we are fitting the inner region of the 
original data, i.e. we did not subtract the outer exponential disc. The 
result is a reasonably good fit, suggesting that the inner r <2 arcsec 
region - where, as noted, the isophotes are quite round - is a sepa- 
rate component (best-fit Sersic parameters: n = 2.35, Rg = 1.60 arc- 
sec, /jp = 14.53). 



2.3.3 Kinematic structure of the photometric bulge 



In Fig. |5]we use the long-slit kinematic data of Heraude au et aE] 
( ,1999) to show an estimate of the local ratio of ordered to random 
stellar motions as a function of radius: Vjp, which is the observed 
stellar velocity deprojected to its in-plane value (assuming an ax- 
isymmetric velocity field), divided by the observed velocity dis- 
persion at the same radius. Even though the data are all inside the 
photometric bulge (r < 50 arcsec), the ratio of Vdp/cr rises above 1 
over much of this region. This is certainly higher than one would 
expect for a classical (kinematically hot) bulge, in which stellar mo- 
tions are dominated by velocity dispersion. An unpublished spec- 
trum with higher S/N from the Hobby-Eberly Telescope (M. Fabri- 
cius, private communication) shows even larger values of Vdp/o" 
for r > 30 arcsec, as well as Vdp/cr > 1 on both sides of the cen- 
tre at r ~ 5-9 arcsec. Our tentative conclusion is that most of the 
photometric bulge is thus a discy, kinematically cool pseudobulge 
(in effect, an inward extension of the disc), similar to that found in 
NGC 3945. We note that the photometric bulge of NGC 3368 has 
also been classified as a pseudobulge by |Drory & Fisherj ( |2007[ l, 
based on morphological features in HST images. Dynamical mod- 
elling of the high-resolution SINFONI data shows that the centre 
of the photometric bulge harbours a kinematically hot component 
(see §4.4), thus confirming the presence of a small classical bulge 
component. 
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Figure 6. Isophotal maps of NGC 3489 on different scales, all with logarithmic intensity scaling. Each successive panel is a zoom of the region outlined with a 
box in the preceding panel, (a) SDSS r-band isophotes, showing the outer disc, (b) //-band isophotes from our WHT-INGRID image, (c) Same, now showing 
the bar oriented almost vertically, (d) Now showing dust-corrected F814W isophotes from the PC chip of the WFPC2 image, (e) Same as previous, showing 
the inner part of the pseudobulge region, (f) Same as previous, but now showing the classical bulge region and the possible nuclear disc. The dashed blue 
square shows the approximate field of view of our 25mas SINFONl observations. 



2.3.4 Orientation and inclination of the galaxy 

Our ellipse fits of the merged SDSS r-band image (black points in 
Fig. [2]l shows a consistent position angle of » 172° for the out- 
ermost isophotes. These ellipse fits extend well outside the star- 
forming outer ring {r ~ 180 arcsec) and are thus unlikely to be 
affected by any intrinsic noncircularity of the ring itself. This po- 
sition angle agrees very well with the kinematic position angles 
determined from both Hi observations (PA x 170°, based on the 
data of |Schneider|1989| as reported by |Sakamoto et al.|1999 i and 
from the Fabry-Perot Ha-l-[Nii] velocity field of Si l'chenko et ak] 
( |2003| see also |Moiseev, Valdes & Chavushyan,2004t . Sil'chenko 
et al. also find a kinematic position angle of 170-175° in the stellar 
kinematics of the inner 2-5 arcsec, from their IFU data. 

The ellipticity of the outer r-band isophotes is « 0.37, with a 
range of 0.34-0.39. A lower limit on the inclination is thus 51°, for 
a razor-thin disc; thicker discs imply higher inclinations. For an in- 
trinsic thickness of c/a = 0.2-0.25, the inclination is / x 53°. This 
is close to the inclination of 50° estimated by [Barbera, Athanas^ 
soula & Garcia-G6mez| pb04[ >, based on Fourier analysis of lFrei| 
et ak| ( |1996) images (note that these images do not extend beyond 



the outer-ring region, and so they might in principle be biased if the 
ring is noncircular). 

An additional, independent estimate of the inclination can be 
had by inverting the Tully-Fisher relation: since we know the ob- 
served H I velocity width and the distance to NGC 3368, we can de- 
termine the inclination needed to make the galaxy follow the Tully- 
Fisher relation. We use the recently published 2MASS Tully-Fisher 
relation of [Masters, Springob & Huchra| pOOS) . For a 2MASS 
AT-band "total" magnitude (Jarrett et al. 2003) of 6.31 (including 



a slight reddening correction from Schlegel, Finkbeiner & Davis| 
|1998[ as given by NED) and a distance of 10.4 Mpc, the absolute 
magnitude is Mi; = -23.795. Using the "Sb" T-F relation from 
table 3 of [Masters et al.| ( |2008[ ), this corresponds to a corrected, 
edge-on velocity width of W^orr = 425 km s"'. For the observed 
width, we use the tabulated value in |Springob et al.| ( |2005| l, which 
is Wc = 324 km s"' ; after applying the recommended correction for 
turbulent broadening (6.5 km s"'), this becomes Wcon- sin i = 317.5 
km s"', and thus i = 48°. 

Taken all this into considerations, we can argue that 
NGC 3368 has a line of nodes with PA » 172° and an inclination 
somewhere between 48° and 55°, most likely a 53°. 
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Figure 7. Isophotal ellipse fits for NGC 3489, plotted on logarithmic (left panel) and linear (right panel) scales. Black points are for the SDSS r-band image, 
while red points are for the INGRID //-band image and purple points are for the dust-corrected WFPC2 F814W image. Major morphological features are 
indicated. Note that due to projection effects, the bar shows up as both a strong twist in the position angle and a minimum in the ellipticity. The dashed line in 
the top panels indicates our adopted position angle (71°) for the galaxy disc. 



2.4 NGC 3489 

2.4. 1 Morphological overview and evidence for a composite 
pseudobulge 

NGC 3489 is structurally somewhat simpler than NGC 3368, with 
only one bar instead of two. The bar itself is rather weak and dif- 
ficult to recognise, because it lies almost along the minor axis of 
the galaxy. Projection effects thus foreshorten it so that it is visible 
primarily due to the abrupt isophote twists, manifesting in the el- 
lipse fits as an extremum in the fitted position angle and a minimum 
in the ellipticity (Fig. |6] and Fig.^. Further outside, the isophotes 
become maximally elongated at r ~ 50 arcsec, and then converge 
to a mean ellipticity of « 0.41 at larger radii. As shown by |Erwin] 
|& Sparke (20031, the ellipticity peak at r ~ 50 arcsec is due to an 
outer ring; the lower ellipticity outside is thus the best representa- 
tion of the outer disc. 



2.4.2 Bulge-disc decomposition 

We do find some evidence that the inner structure of NGC 3489 
is similar to that of NGC 3368 (except for the absence of an inner 
bar in NGC 3489). We start with the global bulge-disc decomposi- 
tion (Fig.[8]l, which results in a Sersic component (the photometric 
bulge) dominating the light at r < 10 arcsec. The isophotes in this 
region are still fairly elliptical - e.g., ellipticity « 0.33 at r ~ 6- 
8 arcsec, which suggests that the photometric bulge is a flattened 
structure. 

As in the case of NGC 3368, we find that the profile of the 
inner photometric bulge can be decomposed into an exponential 
plus a smaller, additional Sersic component (Fig. |9]l. This Sersic 
component dominates the light at r < 2 arcsec. Note, however, that 
the isophotes do become quite elliptical at r ~ 1 arcsec, so we do 
not have as clean a case as in NGC 3368 for a rounder spheroidal 
component. The fit in Fig. |9] shows evidence for a possible nuclear 
excess at r < 0.5 arcsec. This might be evidence for a separate 
nuclear star cluster, similar to that seen in the profile of NGC 2787 
( |Erwin et al.|2003| >. 
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Figure 8. Global bulge-disc decomposition of NGC 3489. The data points 
(circles) are the major-axis cut, combining our ground-based //-band image 
(r > 1.7 arcsec) with the dust-corrected WFPC2 F814W image (scaled to 
//-band). Also shown is the best Sersic + exponential fit to the data and the 
residuals (bottom panel). The Sersic component represents the "photometric 
bulge," which dominates the light at r < 10 arcsec. (As with NGC 3368, we 
argue that most of this is really part of the disc; see text for details.) 



2.4.3 Kinematic structure of the photometric bulge 

There is in addition kinematic evidence that the exponential part of 
the photometric bulge region is kinematically cool and thus a pseu- 
dobulge. In Fig. [To] we plot the local ratio of (deprojected) stellar 
rotation velocity to velocity dispersion. These values are based on 
synthesised long-slit profiles derived from the SAURON and OA- 
SIS velocity and velocity dispersion fields ( [Emsellem et al.|2004[ 
[McDermid et al.||2006| , using slits at PA = 71°. The Vdp/cr ra- 
tio rises to values > 1 at r > 7 arcsec, still within the photometric 
bulge-dominated region, which suggests that the photometric bulge 
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Figure 9. Bulge-disc decomposition of the inner pliotometric bulge of 
NGC 3489. The data points (circles) are the //-band profile from Fig. [s] 
for the inner r < 13 arcsec (with data at r < 1.7 arcsec coming from the 
dust-corrected WFPC2 F814W image). Also shown is the best Sersic + ex- 
ponential fit to the data at r < 13 arcsec and the residuals (bottom panel), 
with the Sersic component representing the classical bulge. 
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Figure 10. Local estimates of the ratio of ordered motion to random motion 
in stellar kinematics along the major axis of NGC 3489, using major-axis 
profiles extracted from the velocity and velocity dispersion fields of the 
SAURON (black) and OASIS (green) observations; the vertical green dot- 
ted lines indicate the seeing of the OASIS observations (FWHM = 0.69 arc- 
sec). Stellar velocities are deprojected to their in-plane values, then divided 
by the observed velocity dispersion values. Ratios < 1 indicate a domi- 
nance of velocity dispersion over bulk rotation. At r < 10 arcsec (vertical 
dashed line), the photometric bulge dominates the light (see Fig. [8}. Since 
Vdp/o" > 1 within this region, the outer part of the photometiic bulge is still 
dominated by rotation, making it a kinematic pseudobulge. 



is at least partly dominated by rotation. As in NGC 3368 Vdp/o" 
drops in the centre but only at smaller radii. Together with the sig- 
nificant flattening in the centre this suggests that the very central 
region of NGC 3489 is structurally different (dynamically colder), 
which is also supported by the dynamical analysis (see §5.5). 



2.4.4 Orientation and inclination of the galaxy 

Lacking the extensive large-scale gas kinematic information that 
was available for NGC 3368, we rely on the isophotes of the outer 
disc to determine the global orientation of NGC 3489. Fortunately, 
apart from the local maximum in ellipticity at r ~ 50 arcsec due to 
the outer ring (see above), the outer disc is fairly well defined, with 
position angle = 71° and a mean ellipticity of 0.41, corresponding 
to an inclination of 55°. 



3 SPECTROSCOPY 

3.1 Data & Data Reduction 

NGC 3368 and NGC 3489 were observed between March 22 and 
24, 2007, as part of guaranteed time observations with SINFONI 
( [Eisenhauer et al."2Q03 Bonnet et al. 2004), an adaptive-optics 
assisted integral-field spectrograph at the VLT UT4. We used 
the ^-band grating and the 3x3 arcsec^ field of view (0.05 x 
0.1 arcsec^ spaxeL') for NGC 3368 and the 0.8 x 0.8 arcsec^ field 
of view (0.025 x 0.0 1 25 arcsec^ spaxeL ' ) for NGC 3489. The total 
on-source exposure time was 80 min for NGC 3368 and 120 min 
for NGC 3489, consisting of 10 min exposures taken in series of 
"object-sky-object" cycles, dithered by a few spaxels. 

The laser guide star (LGS) PARSEC ( |Rabien et aLl[2004l 
[Bonaccini et al.| |2002[ > was used for the AO correction of 
NGC 3368, with the tip-tilt sensor closed on the nucleus of 
NGC 3368 {R = 13.58, B- R = 1.86 within a 3 arcsec diame- 
ter aperture). Although the nucleus itself is just bright enough to 
be used as natural guide star, its shape is rather irregular and not 
pointlike in the /?-band due to the large amounts of dust in the 
nuclear regions and the lack of a strong AGN. Therefore a better 
AO correction was expected from using the LGS instead. The am- 
bient conditions were good and stable, with an average seeing of 
a; 0.6 arcsec in the near-IR. The point-spread function was derived 
by taking an exposure of a nearby star with approximately the same 
/?-band magnitude and B — R colour as the nucleus of NGC 3368, 
using the LGS with the PSF star itself as tip-tilt reference star. The 
FWHM of the PSF is x 0.165 arcsec (see left panel of Fig.[TT]( and 
the achieved Strehl ratio is a; 14%. Due to the time gap between the 
observations of the galaxy and the PSF star the measured PSF shape 
could be different from the PSF during the galaxy observations. We 
compared the surface brightness profile of the SINFONI data with 
the surface brightness profile of an HST NICM0S2 F160W im- 
age, convolved with Gaussians of different widths and found that 
the NICMOS surface brightness profile most closely resembles the 
SINFONI profile for a FWHM ^ 0.165 (see right panel of Fig.[TTJ, 
confirming our PSF measurement. Note that the Gaussian fitted to 
the PSF in the left panel of Fig. [TT] was only used to determine a 
nominal spatial resolution and as a reference for comparison with 
the NICMOS surface brightness profile. A single Gaussian does 
not fit the wings of the PSF, however, the discrepency between the 
PSF and the fit is only ~ 3% in integrated flux. For the dynamical 
modelling we do not use this fit, but the observed image of the PSF 
star. 

NGC 3489 was observed using its nucleus with R - 13.22 
(3 arcsec diameter aperture) as natural guide star for the AO cor- 
rection. A PSF star with a similar magnitude and B — R colour 
was observed regularly in order to determine the spatial resolu- 
tion. The ambient conditions were excellent and stable with a see- 
ing around 0.5 arcsec in the near-IR, resulting in a FWHM of the 
PSF of ~ 0.08 arcsec and a Strehl ratio of 43% (see Fig.[T2l. 
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Figure 11. Left panel: SINFONI PSF derived by observing a star of the same magnitude and colour as the nucleus of NGC 3368 (red dashed line). A Gaussian 
fit is overplotted for comparison (black solid line). Its FWHM is ~ 0.165 arcsec. Right: Comparison of the SINFONI X^-band surface brightness profile with 
the surface brightness profile of an HST N1CMOS2 F160W image convolved with a Gaussian such that the spatial resolution is 0.165". The SINFONI profile 
is shifted such that it matches the NICMOS profile. 




NGC 3489 PSF 
Gaussian fit 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 
r (arcsec) 

Figure 12. PSF of NGC 3489 with a FWHM of ~ 0.08 arcsec. As in Fig. 
|1 l| a Gaussian fit is overplotted for comparison. 



The data reduction was done using the SINFONI data reduc- 
tion package SP RED ijSchreiber et alll[2004{ [Abuter et ar][2006) 
as explained in |Nowak et al. j2008^ . The reduction of the tel- 
luric standard and the PSF reference star was done with the ESO 
pipeline. For the flux calibration we used the telluric standard stars 
Hip 046438 and Hip 085393 with 2MASS Ks magnitudes of 7.373 
and 6.175 respectively as a reference. Fig. |13| shows the flux- 
calibrated images of the two galaxies, collapsed along the wave- 
length direction. 



3.2 Stellar Kinematics in NGC 3368 

The SINFONI data of NGC 3368 were binned using a binning 
scheme with five angular and ten radial bins per quadrant, adopt- 
ing a major-axis position angle of 172°. As in |Nowak et al.| j2008^ 
we used the maximum penalised likelihood (MPL) technique of 
IGebhardt et al.| l [2000a^ to extract the stellar kinematics from the 
first two CO bandheads '2cO(2~0) and '-C0(3-l), i.e. the spec- 
tral range between 2.279 |i.m and 2.340 |i.m rest frame wavelength. 
With the MPL method, non-parametric line-of-sight velocity dis- 
tributions (LOSVDs) are obtained by convolving an initial binned 
LOS VD with a linear combination of template spectra. The residual 
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Figure 13. SINFONI images with overplotted isophotes of NGC 3368 (up- 
per panel) and NGC 3489 (lower panel). The contour levels increase lin- 
early. 
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differences between the resulting model spectrum and the observed 
galaxy spectrum are then calculated. Then the velocity profile and 
the template weights are successively adjusted in order to optimise 
the fit by minimizing the function x\ = + °''^> where a is the 
smoothing parameter and f is the penalty function. The S/N of the 
birmed spectra ranges between 80 and 120 with a mean value of 
~ 1 10. In order to determine the optimal smoothing parameters we 
performed simulations in the same way as in Appendix B of ]Nowak| 
|et al . ( 2008 1, but tailored to our dataset. For a galaxy with a velocity 
dispersion around 100 km s"', a velocity bin width of ~ 35 km s"' 
and the mentioned S/N a smoothing parameter a x 5 is appro- 
priate. As kinematic template stars we chose four K and M giants 
which have about the same intrinsic CO equivalent width (EW) as 
the galaxy (12-14 A, using the EW definition and velocity disper- 
sion correction from |Silge & Gebhardt|2003| l. The uncertainties on 
the LOSVDs are estimated using Monte Carlo simulations ( |Geb-| 
|hardt et ar]|2000al l. First, a reference galaxy spectrum is created 
by convolving the template spectrum with the measured LOSVD. 
Then 100 realisations of that initial galaxy spectrum are created by 
adding appropriate Gaussian noise. The LOSVDs of each realisa- 
tion are determined and used to specify the confidence intervals. 
We verified that the error bars are correct by checking that the S/N 
in the simulated spectra corresponds to the S/N measured from the 
galaxy spectrum. As shown in |Nowak et al.|p008[ >, possible biases 
in the measured LOSVDs are always smaller than the statistical 
errors. 

For illustration purposes we fitted Gauss-Hermite polynomi- 
als to the LOSVDs. Fig. [14] shows the two-dimensional fields of 
V, <r and the higher-order Gauss-Hermite coefficients and h4, 
which quantify the asymmetric and symmetric deviations from a 
Gaussian velocity profile l |Gerhard||1993[ [van der Marel & Franx| 



[1993| >. The major-axis profiles are shown in Fig. [21 



The velocity field of NGC 3368 shows a regular rotation about 
the minor axis. The average, luminosity-weighted cr within the to- 
tal SINFONI field of view is 98.5 km s"'. A central cr-drop of 
7% is present within the inner ~ 1 arcsec, well inside the region 
of the classical bulge component, cr-drops are not uncommon in 
late-type galaxies and are usually associated with nuclear discs or 
star-forming rings (e.g. [Wozniak et al.|2003| [Peletier et al.|2007| 
[Comeron et al.|2008) . These could be formed e.g. as a result of gas 
infall and subsequent star formation, but as no change in ellipticity 
is found in the centre, such a disc would have to be very close to 
face-on. A cr-drop does not imply the absence of a SMBH if the 
centre is dominated by the light of a young and kinematically cold 
stellar population (see also the discussion in §4.4). |Davies et al.| 
( |2007 1 observed cr-drops in a number of strongly active galaxies. 
In these AGN the mass of the central stellar component was ~ 10 
times that of the SMBH, so no outstanding kinematic signature 
would be expected. Another example is the velocity dispersion of 
the Milky Way, which apparently drops in the central 100 pc, and 
only rises in the inner 1-2 pc (see Figure 9 of |Tremain e et al. 2002). 
Finally, a central cr-drop has been found in NGC 1399 ^Gebhardt| 
|et al.|2007[ |Lyubenova et al.|2008| l, where it has been interpreted 
as a signature of tangential anisotropy. 

The velocity dispersion in quadrants 2 and 3 is smaller than in 
quadrants 1 and 4. A possible explanation for that behaviour could 
be the substantial amounts of dust in the central regions (Fig.|15[l, 
although the effect of the dust in the A'-band is relatively weak. 
The HST WFPC2 B- I colour map (Fig.[T5j shows that within the 
SINFONI field of view the dust extinction is largest in quadrants 1 
and 2. Quadrant 4 is moderately affected while quadrant 3 seems to 




Figure 15. HST WFPC2 B- I colour map of NGC 3368. Indicated is the 
major axis as a solid line (PA= 172°) along with the PA= 5° slit orientations 
used by |Heraudeau et al.|{T999) and |Vega Beltran et al.|j2001) . 



be relatively dust-free. We will further discuss the asymmetries in 
§4.5. 

In the near-IR the presence of dust should have a much smaller 
effect on the kinematics than in the optical, therefore the asym- 
metry should be much stronger in the kinematics measured using 
optical absorption lines, if dust is responsible for the asymmetry. 
Longslit kinematics (v and cr) at PA= 5°, measured from opti- 
cal spectra using Fourier-Fittin g or FCQ jBender, Saglia & GeF 



hard|1994), a re available from |Heraudeau et al.H1999^ and jVega 



Beltran et al.| ( |200T| ). Two-dimensional kinematics have been mea 



sured by ISil'chenko et al.| (^2003 ) (see also |Moiseev et al.||2004[ > 
with the Multi-Pupil Field Spectrograph (MPFS) at the Russian 
6 m telescope in the optical using a cross-correlation technique. 
The spatial resolution of the optical data is between 1.4 arcsec and 
3.0 arcsec. The velocities of the different authors are in good agree- 
ment with each other and with the SINFONI velocities consider- 
ing the different seeing values. The optical velocity dispersions are, 
however, significantly larger than those measured with SINFONI. 
They are on average around 130 km s"' for the longslit data and 
~ 150 km s"' for the MPFS data. There are a number of possible 
causes for such a discrepancy. The authors used different correla- 
tion techniques, slightly different wavelength regions and different 
templates. A difference between optical and AT-band cr measure- 
ments was also found by |SiIge & Gebhardt[ l |2003[ > for a sample 
of galaxies and they suggested that this might be caused by strong 
dust extinction in the optical. But weak emission lines could also 
alter the absorption lines and thus the measured kinematics. As in 
the SINFONI data, a velocity dispersion asymmetry is also present 
in all optical datasets, as well as a velocity asymmetry. The velocity 
dispersion of |Moiseev et al.||2004[ l is enhanced in the entire region 
west of the major axis, where also the majority of the dust is located 
(Fig.|15[l. However, when comparing the extinction along the loca- 
tion of the longslits of |Heraudeau et ar] ( |1999| l and |Vega Beltran| 
|et al.|pOOT; ) with the according velocity dispersion, there seems to 
be no correlation. Thus it is not clear whether and in what way dust 
influences the velocity dispersion in NGC 3368. 

Another explanation for the asymmetry could be lopsided- 
ness, which is common in late-type galaxies. Possible mechanisms 
which could cause lopsidedness include minor mergers, tidal in- 
teractions and asymmetric accretion of intergalactic gas ( [Boumaud[ 
|et al.|2003] >. As the large-scale stellar and gas velocity fields and gas 
distributions jSil'chenko et al.|20 03 ; Haan et al. 2008 1 are rather 
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Figure 14. Two-dimensional stellar kinematics (v, tr, and /14) of NGC 3368. 
velocity map (upper left). 

regular, a recent merger or collision with another galaxy seems un- 
likely. Accretion of gas from the intergalactic HI cloud is a more 
likely scenario ( |Schneider|1989||Sil'chenko et al.|2003l > and could 
be a possible explanation for the presence of molecular hydrogen 
clouds close to the centre (see below). However, there seems to be 
no lopsidedness in the A'-band photometry, as any distortions of 
the isophotes can plausibly be explained by dust. The molecular 
gas distribution on the other hand is very disordered in the central 
~ 200 pc (see below and |Haan et al.|2009 l|. Thus if the gas mass 
differences between different regions of the galaxy would be large 
enough, they could be a plausible explanation for the distorted stel- 
lar kinematics. However, as shown in ^4.5[ the molecular gas mass 
is small compared to the dynamical mass and is thus unlikely to 
have a significant eff^ect on the stellar kinematics. 

Central lopsidedness like an M31-like nucleus or otherwise 
off-centred nuclear disc ( Bender et al.|2005{|Jog & Combes|2009^ 
could, if the resolution is just not high enough to resolve the disc as 
such, leave certain kinematical signatures like a slightly off-centred 
(T-peak or -drop. On the other hand we see velocity asymmetries 
out to r ~ 20 arcsec, which is way too large to be explained by an 
M31-like nuclear disc. 

In principle, the outer and inner bars could cause asymmetries 
in the stellar kinematics. However, the SINFONI field of view is 
located well inside the inner bar, and the only changes in velocity 
dispersion associated with inner bars which have been observed 
are symmetric and take place at the outer ends of inner bars fde' 
[Lorenzo-Caceres et al.|20()8t . 




Major axis, minor axis and the numbering of the quadrants are indicated in the 



3.3 Gas Kinematics in NGC 3368 

In NGC 3368 the only emission lines detected arise from molecular 
hydrogen H2. The strongest line is l-OS(l) at A = 2.1218 \im. To 
determine the flux distribution and velocity of the H2 gas we fitted 
a Gaussian convolved with a spectrally unresolved template pro- 
file (arc line) to the continuum-subtracted spectrum ( ,Davies et al.| 
|2007| l. The parameters of the Gaussian are adjusted such that they 
best fit the data. Fig.[T6]shows the flux distribution and the velocity 
field of H2 l-OS(l). As the S/N of the H2 emission is very low in 
some regions we binned the data using adaptive Voronoi binning 
( [Cappellari & Copin||2003| to ensure an approximately constant 
S/N and thus a robust velocity measurement in each bin. This bin- 
ning scheme is different from the radial and angular binning used 
to measure the stellar kinematics (see Fig. |14[ l, which is appropri- 
ate for the dynamical modelling procedure. The flux distribution 
of the gas is different and more complex than the flux distribution 
of the stellar light, thus it would be inappropriate to use the same 
binning as for the stars for the purpose of S/N adjustment of the 
gas emission. The most striking feature seen in Fig.[T6]are the two 
clouds of H2 gas, located ~ 0.36" and ~ 0.72" north of the photo- 
metric centre. These two clouds are kinematically decoupled from 
the remaining H2 gas distribution and seem to move in opposite di- 
rections. Their projected sizes are approximately 25 pc and 20 pc 
FWHM. The H2 distribution outside these two clouds is relatively 
smooth. Its kinematic position angle, measured using the method 
described in Appendix C of |Krajnovic et al.]j2006f , is ~ 171° and 
thus agrees with the stellar kinematic position angle. The gas ve- 
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Table 3. Mean near-IR line strength indices in A of NGC 3368 (3x3 arcsec 
apeilure) and NGC 3489 (0.8 X 0.8 arcsec aperture). The corresponding rms 
is given in brackets. 





NGC 3368 


NGC 3489 


Nal 


4.30(0.16) 


3.30 (0.36) 


Cal 


2.40 (0.29) 


1.53 (0.50) 


FelA 


1.55 (0.07) 


1.25 (0.23) 


FelB 


0.91 (0.10) 


0.76 (0.21) 


CO 


17.90 (0.27) 


17.70 (0.72) 



locity follows the rotation of the stars within a radius of ~ 0.5"; 
outside that radius it rotates faster, reaching rotation velocities up 
to - 100 km s"'. 

The origin of the H2 clouds is unclear. H2 emission line ra- 
tios can help to distinguish between different excitation mecha- 
nisms like shock heating, X-ray illumination or UV fluorescence 
( |Rodn'guez-Ard ila, Rifi^el & Pastoriza'2005 1. Table[2]shows the Hj 
line ratios for difi'erent regions (the total SINFONI field of view and 
the two H2 clouds). They indicate that the H2 gas is thermalised in 
all regions of the SINFONI field of view (cf. figure 5 in |Rodrfguez-| 
lArdila et al.|2005l l. 

The ratio 2-lS(l) 2.247 ^m/l-OS(l) 2.122 ^.m and therefore 
the vibrational excitation temperature jReunanen, Kotilainen & Pri^ 
|eto|2002] l is smaller in the H2 clouds than in the regions where H2 
is evenly distributed. 



3.4 Line Strength Indices for NGC 3368 

The stellar populations of NGC 3368 have been analysed by 
ISil'chenko etaTlpOOS) and |Sarzi et allpOOS) using optical spectra. 
Both found that a relatively young stellar population with a mean 
age of around 3 Gyr dominates the central region. Towards larger 
radiilSil'chenko et al.H2003|l found a strong increase in age. 



Fig. 17 shows the near-IR line indices Nal and CO measured 



in the same way as in 



Nowak et al. 



1 2008b, using the definitions 



of |Silva, Kuntschner & Lyubenova ^0081. The average values in 



side the SINFONI field of view are listed in Table |3] They differ 
significantly from the relations between Nal or CO and cr found 
by |Silva et al.| ( [2008 1 for early-type galaxies in the Fornax clus- 
ter which may be, as in the case of Fornax A, probably due to the 
relatively young age of the stellar population. Younger populations 
seem to have larger Nal at equal a than old stellar populations in 
the galaxy samples of |Silva et ar] ( |2008| l and |Cesetti et al.|p009l l. 
However, no such trend is obvious for the CO index, so the dif- 
ference seen here could be due to other aspects like metallicity or 
galaxy formation history. The measured average indices have val- 
ues which are quite similar to those found in the centre of Fornax A 
( [Nowak et al.|2008] l, which could indicate that the stellar popula- 
tions are quite similar in terms of age and metallicity. However, 
the interpretation of Na I, Ca I and <Fe I> must be done bearing 
in mind that these features always include significant contributions 
from other elements ( |Silva et al.|2008[ l. 

The radial distribution of the line indices (Fig. [TtJ shows a 
slight asymmetry, similar to the kinematics. The two quadrants with 
the smaller cr have larger CO EWs and smaller Na I EWs than the 
other two quadrants. In addition there seems to be a strong negative 
gradient in Nal and a moderately strong negative gradient in CO. 
Cal and <FeI> are approximately constant with radius. A small 



central drop is present in most indices, which could indicate the 
presence of weak nuclear activity (see [Davies et al.|2007^ . 



3.5 Stellar Kinematics in NGC 3489 

The NGC 3489 data were binned in the same way as the NGC 3368 
data, with identical angular bins and nine somewhat smaller radial 
bins due to the higher spatial resolution. A position angle of 71° 
was used. The mean S/N is 70, and a smoothing parameter of a = 
8 was used. As kinematic template stars we chose four K and M 
giants with an intrinsic CO equivalent width in the range 13-15 A. 

Fig. [18] shows the v, cr, and /14 maps of NGC 3489. The 
kinematics is similar to that of NGC 3368 in some aspects. It 
is clearly rotating about the minor axis, though stronger than 
NGC 3368. The velocity dispersion also drops towards the centre 
by around 4%, but then has a tiny peak in the central bins. The av- 
erage cr in the total 0.8 x 0.8 arcsec^ field of view is 91 km s"' . The 
values clearly anticorrelate with v. /14 is on average small and 
negative. Along the major axis it is positive in the outer bins and 
negative in the inner bins. No asymmetry is present in cr, but the 
velocity is slightly asymmetric. It increases strongly on the reced- 
ing side and then remains approximately constant at r > 0.05 arc- 
sec, whereas on the approaching side the slope is less steep and 
an approximately constant velocity is reached much further out 
at r > 0.2 arcsec. The major-axis profiles are shown in Fig. [26] 
Note that the central velocity bin is omitted in that plot. Despite the 
presence of strong dust features in optical images, the kinematics 
is in comparatively good agreement with the 2D SAURON ( |Em-[ 
[sellem etal. 2004 1 and OASIS ( McDermid et al.|2006 ) kinematics, 
though due to the high spatial resolution and the very small field 
of view of our data a direct comparison with seeing-limited data is 
not easy. The SINFONI velocity field seems to be fully consistent 
with the optical velocity. The average SINFONI velocity dispersion 
is smaller than the central SAURON and OASIS cr. The central 
SAURON and OASIS is significantly larger than the SINFONI 
values, and the anticorrelation of and v seems to be less strong 
in general and essentially non-existent in the central arcsecond of 
the OASIS data. 

Longslit kinematics is available from [Caon, Macchetto & Pas^ 
[toriza|pb00^ , who used a cross-correlation technique to determine 
V and cr. Their cr is much larger (117 km s"' in the central pix- 
els of the major-axis longslit), which could be due to their cross- 
correlation technique or template mismatch. In addition the cr- 
profile seems to be slightly asymmetric, but only at large radii, 
where the errors are large. 



3.6 Line Strength Indices for NGC 3489 

The stellar populations of NGC 3489 in the centre have been anal- 
ysed by [Sarzi et al. ( 2005 ^ using optical HST STIS longslit spectra 
and by [McDermid et al.[ ^2006^ using OASIS integral-field data. 
[Sarzi et al.[j2005> obtained a mean age of about 3.1 Gyr in the cen- 
tral 0.2 X 0.25 arcsec by fitting stellar population synthesis models 
to the spectra, assuming solar metallicity. [McDermid et al.[ p006[ > 
obtained a mean age of 1.7 Gyr in the central 8x10 arcsec'^ with an 
age gradient down to ~ 1 Gyr towards the centre from the analysis 
of Lick indices. These two values are more or less in agreement 
when taking into account the measurement errors and that the [Sarzi[ 
[et al.[ ( [2005) value would decrease when considering the metallicity 
increase to supersolar values in the centre measured by |McDermid[ 
[et al.j (|2006jl. Another possibility is that the central ~ 0.2 arcsec. 
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Figure 16. Left panel: H2 (1-0S(2)) emission in the centre of NGC 3368. Right panel: H2 gas velocity. The dashed line indicates the major axis. The small 
inset shows the pseudo-longslit gas velocity profile along the major axis (black points) in comparison with the major-axis stellar velocity (red points). The 
isophotes of the stellar emission are overlaid in both panels. 



Table 2. H2 l-OS(l) 2.12iam emission line fluxes and H2 line ratios. 



Region 


l-0S(l)(10-'5ergs-' cm-) 


1-0S(3) 
l-OS(l) 


1-0S(2) 
l-OS(l) 


l-OS(O) 
l-OS(l) 


2-lS(l) 
l-OS(l) 


Tv.b (K) 




Mf[°l'' (10^ Moj^ 


cloud 1 


1.66 


1.35 


0.56 


0.20 


0.12 


2280 


9.13 


2.23 — 0.21 


cloud2 


0.83 


1.30 


0.51 


0.21 


0.11 


2203 


4.57 


1.12 — 0.10 


total 


21.40 


1.69 


0.98 


0.26 


0.17 


2740 


117.77 


28.72 — 2.65 



The first value gives the cold H2 mass as estimated from hot H2 using the conversion of |Mueller Sanchez et al. J (2006). 
The second value gives the real cold gas mass after calibration against direct mass measurements from CO (J = 1 - 0) at 
r = 1 .5 arcsec jSakamotoetal.|1999| . 



(a) NGC 3368 (b) NGC 3489 




r (arcsec) r (arcsec) r (arcsec) r (arcsec) 

Figure 17. Near-IR line indices Nal and CO in A for (a) NGC 3368 and (b) NGC 3489. The upper row shows the two-dimensional index fields, with the 
position-index diagrams in the bottom row. Individual values are plotted in grey, while overplotted in red are the mean indices and their rms for bins belonging 
to semicircles in the receding half of the galaxy (negative radius r) or the approaching half (positive radius r). 
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NGC 3489 
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Figure 18. Two-dimensional stellar kinematics (v, tr, and /14) of NGC 3489. Major axis, minor axis and the numbering of the quadrants are indicated in the 
velocity map (upper left). 



which are unresolved by |McDermid et ar] ( |2006[ l, contain an older 
stellar population. 

The near-IR absorption-line indices Nal and CO are shown in 
Fig.|17|3. They seem to be, like the stellar kinematics, axisymmet- 
ric. As in the case of NGC 3368 there is a clear negative gradient 
in both indices, which could mean an age or a metallicity gradi- 
ent, or a combination of both. The other indices, Cal and <FeI>, 
are largely constant. The average line strength indices within the 
0.8 X 0.8 arcsec^ field of view are given in Table [5] The average 
CO line strength is very similar to the value found in NGC 3368. 
All other measured indices are slightly smaller than in NGC 3368. 
This seems to be generally in agreement with the results of ,Sarzi, 
|et al.|p005) , who found similar mean ages and populations in both 
galaxies. A small central drop is present only in Nal, implying that 
nuclear activity must be extremely weak or absent. 



4 DYNAMICAL MODELLING OF NGC 3368 

For the dynamical modelling we make use of the |Schwarzschild| 
( |1979^ orbit superposition technique: first the gravitational poten- 
tial of the galaxy is calculated from the stellar luminosity density 
V and trial values for the black-hole mass M. and the mass-to-light 
ratio T. Then an orbit library is generated for this potential and a 
weighted orbit superposition is constructed such that it matches the 
observational constraints. Finally everything is repeated for other 
potentials until the appropriate parameter space in M. and Y is sys- 
tematically sampled. The best-fitting parameters then follow from a 



^^-analysis. The deprojected luminosity density is a boundary con- 
dition and thus is exactly reproduced, while the LOSVDs are fitted. 



4.1 Construction of the Stellar Luminosity Profile 

For dynamical modelling purposes, we need an appropriate 
surface-brightness profile and an appropriate ellipticity profile, 
along with an assumption of axisymmetry. While simply using the 
results of ellipse-fitting may be valid for an elliptical galaxy, where 
the approximation that the galaxy is a set of nested, axisymmetric 
ellipsoids with variable axis ratio but the same position angle is of- 
ten valid, a system like NGC 3368, with two bars, dust lanes, and 
spiral arms, is clearly more complicated. Such a complex structure 
also makes it important to allow a Y gradient in order to account for 
stellar population changes. This can be conveniently approximated 
by using more than one component, where each component has its 
own Y. 

We model the luminosity distribution of NGC 3368 as the 
combination of two axisymmetric components: a disc with fixed 
(observed) ellipticity = 0.37, which by design includes both inner 
and outer bars and the discy pseudobulge; and a central "classical" 
bulge of variable (but low) ellipticity. Thus, we assume that the bars 
can, to first order, be azimuthally "averaged away." 

The surface brightness profile of the disc component is not as- 
sumed to be a simple exponential. Instead, it is the observed surface 
brightness profile of the entire galaxy outside the classical bulge, 
out to r = 130 arcsec, along with an inward extrapolation to r = 0. 
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We base this profile on ellipse fits with fixed ellipticity and posi- 
tion angle (e = 0.37, PA = 172°) to the A'-band image of |Knapenl 
|et al.|p003] l, with the inner r < 3.7 arcsec based on the exponential 
component of our inner bulge-disc decomposition (Fig.|4|. (Com- 
parison of profiles from the dust-corrected NICM0S2 image and 
the Knapen et al. image shows that seeing affects the latter only for 
r < 2 arcsec, which is already within the region where the classi- 
cal bulge affects the profile.) Inspection of both this profile and a 
similar fixed-ellipse profile from the NICMOS2 image shows that 
the classical bulge begins to affect the profile only for r < 3.7 arc- 
sec. Consequently, the disc profile for r < 3.7 arcsec is the in- 
ward extrapolation of the exponential component from our inner 
Sersic+exponential decomposition (fig = 12.75, h = 5.28 arcsec; 
Fig.g. 

To generate the profile of the classical bulge, we assume, fol- 
lowing the iimer decomposition discussed above, that the light in 
the inner r < 8 arcsec is the combination of an inner exponen- 
tial and the classical bulge (Fig. |4j. We generated a model image 
with the same size as the NICMOS2, containing a 2D exponential 
model for the inner disc, which we subtracted from the NICM0S2 
image. The residual image is assumed to contain light from the 
classical bulge only; we then fit ellipses to this image. This allows 
for possible variations in the classical bulge's ellipticity and, per- 
haps more importantly, uses the observed surface brightness profile 
at the smallest radii, rather than an analytic fit. Finally, we generate 
an extension of this bulge profile out to the same outermost radius 
as the disc profile (i.e., well outside the NICM0S2 image) by fit- 
ting a Sersic function to the classical-bulge profile, and assuming 
a constant ellipticity of and the same PA as for the outer disc at 
large radii. 

The program of |Magorrian| ( | 1 999) was used for the deprojec- 
tion assuming that all components are axisymmetric. Both com- 
ponents, the disc and the classical bulge, were deprojected for an 
inclination / = 53° as obtained from the photometry (see §2.3.4), 
and for a few nearby values between 52° and 55°. No shape penalty 
was applied. The simplest assumption for the form of the stellar 
mass density p, is then p, = Tbuigc • V(,i,igc + Tjisc • vjisc, where v is the 
deprojected luminosity density and Tbuigc and Tdisc two constants to 
be determined jPavies et al. [2006) 1. The assumption of a constant 
T is approximately true for the central part of the galaxy where we 
have kinematic data (dark matter does not play a significant role). 

As a further test, we deprojected the two components from 
the global bulge-disc decomposition (i.e. the outer exponential disc 
and a Sersic fit to the photometric bulge region, as shown in Fig. 
[3](. The resulting shape of the luminosity profile is very similar to 
the profile obtained from the inner bulge-disc decomposition of the 
photometric bulge. The global profile is offset to smaller luminosi- 
ties, as the global decomposition does not fully account for the light 
in the classical bulge component. This means that the SMBH mass 
estimates would be larger for models based on the global bulge- 
disc decomposition (for constant T) compared to the mass estimate 
derived from models using the inner bulge-disc decomposition. 



4.2 Dynamical Models 



As in |Nowak et aT] (|2007[ |2008^ we use an axisymmetric code 
( |Richstone & Tremaine|1988||Gebhardt et al.|20 00a 2003: |Thoma"sl 
[eraL][200? l to determine the mass of the SMBH in NGC 3368. 
This method has been successfully tested on the maser galaxy 
NGC 4258 in |Siopis et jJ^lpOOg) , who obtained the same mass for 
the black hole as determined from maser emission. 

We allow for different mass-to-light ratios in the classi- 



cal bulge region and the region further out, respectively. Radial 
changes in the mass-to-light ratio can bias the derived BH mass 
if not taken into account properly (e.g., |Gebhardt & Thomas|2009| l. 

Using an axisymmetric code for a barred and therefore obvi- 
ously non-axisymmetric galaxy might be debatable. Tho mas et al.| 
(|2007j have found that axisymmetric dynamical models of ex- 
tremely triaxial/prolate systems are in danger of underestimating 
the luminous mass in the centre. Since BH mass and central lumi- 
nous mass are partly degenerate, this could result in an overesti- 
mate of the BH mass. For two-component models the situation is 
even more complex. If the triaxiality only affects the outer Tdisc (as 
in the case of NGC 3368), then a corresponding underestimation 
of the outer T^is;. could translate into an overestimation of the inner 
Tbuigc, which in turn would imply a bias towards low BH masses. 
Detailed numerical simulations of barred galaxies are required to 
investigate such possible biases. 

However, in this case the axisymmetric models can be justified 
as we only model the central part of the galaxy. Near the SMBH the 
potential is intrinsically spherical and strong non-axisymmetries 
are unlikely. Also, there is little evidence that non-axisymmetric 
bar orbits dominate the observed region in projection (see §4.5). 

We use only the SINFONI data for the modelling. The four 
quadrants are modelled separately in order to assess the influence 
of deviations from axisymmetry. These four independent measure- 
ments of the SMBH mass should agree within the observational 
errors, if the galaxy is axisymmetric. If not, then the systematic 
differences from quadrant to quadrant provide an estimate for the 
systematic errors introduced by assuming axial symmetry. 

Since the observed ellipticity of the classical bulge is affected 
by strong dust lanes, the ellipticity is slightly uncertain (see §2.3.1). 
We ran dynamical models for two different deprojections, one ob- 
tained for a bulge ellipticity 6 = 0.0 and the other for 6 = 0.1. The 
models yield the same mass-to-light ratios and black hole masses 
and in the following we only discuss the case 6 = 0.0. 

In order to find out whether the results depend on the assumed 
inclination of the galaxy, we run models for four different inclina- 
tions around the most likely value of 53°. 

We do not apply regularisation to our models, because the ex- 
act amount of regularisation is difficult to determine due to the lack 
of realistic, analytical models of disc galaxies with black holes. 

The SINFONI observations mainly cover the classical bulge 
region, so the disc T can only be weakly constrained. It could be 
better constrained if we included other kinematic data extending 
further out, but this has several disadvantages. The inconsistencies 
between the SINFONI and the optical measurements from the lit- 
erature mean that the models could have difficulties fitting the dif- 
ferent kinematic datasets reasonably well at the same time. In ad- 
dition, the non-axisymmetries due to the two bars would be more 
noticeable at large radii, and a dark halo would become important. 



4.3 Results 
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The results for ; = 53° and 6 = 0.0 are shown in Fig. 
a function of M. and the T of one component, marginalised over 
the other component's T), and for all inclinations in Fig. |20| (total 
X' as a function of one of the three parameters M., Tbuigc, Tjisc, 
marginalised over the other two parameters). The best-fitting values 
with 3cr errors are listed in Table |4] for all four inclinations. Each 
quadrant provides an independent measurement of M. if deviations 
from axisymmetry do not play a big role. This seems to be the case, 
as the resulting best-fitting values for M., Tbuigc and Tjisc agree very 
well wiffiin < 2cr between the four quadrants. The mean black hole 
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mass for the four quadrants (/ = 52 - 55°) is (M.) = 7.5 X 10* Mq 
(rms(M.) = 1.5 X 10'' Mq). Note that the average of all 3 a errors 
given in Table |4] derived from the x' analysis divided by three is 
1.6 X lO*" Mg and thus approximately equal to the rms error. 

The resulting black hole mass does not depend much on the 
particular choice of the mass-to-light ratio of the disc Tjisc, but de- 
creases for increasing Tbuijc. As shown in Fig.|20] the results also do 
not change systematically with the inclination. This shows that the 
inclination cannot be constrained better by dynamical modelling 
than by a thorough analysis of photometric data, and that a very 
precise knowledge of the inclination is not necessary for dynami- 
cal modelling purposes. 



4.4 Evidence for a black hole in NGC 3368 

V, cr, hi and of the best model at / = 53° and e = 0.0 (major axis) 
are shown in Fig. |2T]for all quadrants. The corresponding best fit 
without black hole would be very similar, which is why we choose 
to plot the differences in instead (see below and Fig. 



22 i. The 



similarity between the black-hole model and the one without black 
hole in the lower order velocity moments, as well as the cr-drop 
towards the centre raise the question where the dynamical evidence 
for the black hole comes from. 

; shows the;jf- difference between the best-fitting model 



Fig. 
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without a black hole and the best-fitting model with a black hole 
for all LOSVDs in all four quadrants. The fit with black hole is 
generally better in 132 of the 180 bins. The largest^- differences 
appear in the LOS VD wings, where the model with black hole gives 
less residuals. As a different illustration of what is shown in Fig. 

shows the increase in Ax^ summed over all angular and 
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Fig 
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velocity bins as a function of radius. Thus for each quadrant five 
(angles) times 21 (velocities) Ax~ values were added at each radius. 
The largest Ax~ increase occurs in the central ~ 0.2 - 0.3 arcsec 
(i.e. ~ 4rsoi). At larger radii the Ax^ increase is less strong. Tt^ige 
is large for the best models without black hole, which can worsen 
the fit in the outer data regions. Therefore it is not surprising that 
improvements of the fit appear at all radii. The total Ax^, summed 
over all LOSVDs, between the best-fitting model without a black 
hole and the best-fitting model with a black hole is given in the last 
column of Table |4] 

The asymmetry of the data is reflected in the error bars. For 
quadrants 1 and 4, which are the quadrants with the higher veloc- 
ity dispersion, the error bars are much larger than for quadrants 2 
and 3. In 2 and 3 the no-black hole-solution is excluded by > 5cr, 
whereas for 1 and 4 it is only excluded by ~ 3.6<t. 



4.5 Discussion 

Upper limits for M. in NGC 3368 have been measured from cen- 
tral emission-line widths by |Sarzi et al.| ( (2002l l and |Beifiori et al.| 
(20091, who obtain 2.7 x 10'' Mq (stellar potential included) and 



4.8 X 10^ Mq (stellar potential not included), respectively. Based 
on the M.-cr relation of Tremaine et al. |( |2002| l, a black hole with 
a mass between 8 x 10* and 2.5 x 10^ Mg would have been ex- 
pected, depending on which cr measurement is used. From the 
relation between M. and .K'-band luminosity ( [Marconi & Hunt] 
|2003 l we would have expected a very high black-hole mass of 
9.2x 10' Mq if it correlates with the total (photometric) bulge lumi- 
nosity (Mjf = -23.42), or a very small mass of only 1.5 x 10* Mg 
if it correlates with the classical bulge luminosity (Mx = -19.48). 
The stellar mass within the sphere of influence rsoi = 




0.4 0.6 
r [arcsec] 

Figure 23. x'^ difference between the best-fitting model without black hole 
and the best-fitting model with black hole summed over all angular and 
velocity bins as a function of distance from the black hole for all four quad- 
rants of NGC 3368. 



GM./a^ » 0.07 arcsec is a 5.7 x 10* Mg. If the best-fitting mass 
for the black hole of M. = 7.5 x 10* Mg were entirely composed 
of stars, the mass-to-light ratio (Tbuigc) ~ 0.41 within rsoi would 
increase to 0.95. This would be typical for an older stellar popula- 
tion (~ 4-7 Gyr for a Salpeter IMF and ~ 10-1 1 Gyr for a Kroupa 
IMF at solar metallicity, using the models o f|Maraston|1998||2005l l. 
However, this would strongly conflict with| Sarzi et al.H2005 l), who 
find that a 1 Gyr old population dominates, with some contribu- 
tions from older and younger populations, resulting in a mean age 
of 3 Gyr. 

As mentioned in §4.2, bar orbits crossing the centre could in 
principle produce deviations from axisymmetry. Non-axisymmetric 
structures such as a prolate central structure can be recognised in 
the kinematics as a low-cr, high-/i4 region if seen edge-on, or as a 
high-tr, low-/?4 region if seen face-on ( [Thomas et al.|[2007] l. This 
could also bias the reconstructed masses. Strictly speaking, this is 
only valid for N-body ellipticals with a central prolate structure. 
Simulations of bars by [Bureau & Athanassoula| ^005 ) resulted in 
variable a and /14, depending on the projection of the bar. How- 
ever, these variations were always symmetric with respect to the 
bar. Thus we would expect symmetry between quadrants where the 
bar appears. This is not the case, therefore it is unlikely that bar 
orbits in projection disturb the central kinematics significantly. 

Dust could in principle influence the kinematics and produce 
distortions or asymmetries, though no clear correlation with the 
dust distribution could be found in §3.2. According to [Baes et al.[ 
l!2003 1, however, dust attenuation should not affect moderately in- 
clined galaxies significantly. 

The mass of the hot H2 in the clouds and in the total 
field of view was estimated via Mh, = 5.0875 x 10"D^/i_os(i) 
( [Rodrfguez-Ardila et al.[[2005[ l an d the total cold gas mass vi a 
log(^Li_os(i)A/S',"') = -3.6 ± 0.32 jMueller Sanchez et al.| 



20061, 
T^) 



where D is the distance of the galaxy in Mpc, 1 is tfte J:l2 (1 - U) 
flux and L is the Hi (1 - 0) luminosity. They are listed in Ta- 
ble |2] Note that the latter conversion has large uncertainties and 



M^" spans at least two orders of magnitude 



that the ratio Mil"' 
m 

(~ 10"'-10"^). In addition, M]g' : Mg** depend s on the far-infrared 
colour /^(60^ra)//,,(100^m) ( [Dale et al.|2005) , which may help to 
place tighter constraints on the ratio. With a far-IR colour of 0.35 
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NGC 3368 




M. (10* Ms) M. (10' Mq) M. (10' Mg) M, Mq) 



Figure 19. Models for NGC 3368 calculated for the quadrants 1-4, an inclination (' = 53° and a classical bulge ellipticity e = 0.0. For each quadrant 
A^o ~ ~^min plotted as a function of the black hole mass M, and the ^T-band mass-to-light ratios Tbuisc (top row) and Tjisc (bottom row). The coloured 
regions are the l-Scr confidence intervals for two degrees of freedom. Each model that was calculated is marked as a black dot; the best-fitting model is 
encircled by a white ring. 



Table 4. Resulting black hole masses and K-hand mass-to-light ratios Tbuigc and Tjisc of NGC 3368. The lower and upper 
3o" limits are given in brackets. The total of the best model with black hole and the;j'^ difference between the best model 
without black hole and the best model with black hole are given in the last two columns. 



Inclination 


Quadrant 


M. [10* Mq] 


"^bulge 




^ min 


'^''^noBH-BH 


52° 


1 


8.0 (3.0, 11.0) 


0.45 (0.30, 0.60) 


0.50 (0.20, 0.65) 


137.353 


26.093 




2 


8.0 (4.0, 10.0) 


0.40 (0.30, 0.45) 


0.35 (0.30, 0.50) 


228.051 


32.828 




3 


8.0 (3.0, 10.0) 


0.40 (0.25, 0.50) 


0.45 (0.25, 0.65) 


185.363 


29.471 




4 


9.0(1.0, 13.0) 


0.50 (0.30, 0.55) 


0.30 (0.25, 0.65) 


141.261 


21.207 


53° 


1 


7.0(1.5, 11.0) 


0.40 (0.30, 0.65) 


0.45 (0.20, 0.65) 


134.681 


17.285 




2 


8.0 (3.0, 10.0) 


0.45 (0.25, 0.50) 


0.20 (0.20, 0.40) 


225.221 


26.828 




3 


9.0 (4.0, 10.0) 


0.40 (0.35, 0.45) 


0.35 (0.25, 0.50) 


180.291 


36.955 




4 


6.0(1.5, 14.0) 


0.40 (0.25, 0.65) 


0.50 (0.20, 0.65) 


144.920 


16.538 


54° 


1 


8.0 (2.0, 11.0) 


0.50 (0.35, 0.60) 


0.35 (0.20, 0.50) 


131.651 


25.002 




2 


4.0(1.0, 11.0) 


0.45 (0.30, 0.55) 


0.30 (0.20, 0.55) 


234.792 


15.106 




3 


8.0 (3.0, 11.0) 


0.35 (0.30, 0.45) 


0.40 (0.25, 0.65) 


186.925 


41.063 




4 


6.0(1.0, 14.0) 


0.55 (0.30, 0.70) 


0.30 (0.20, 0.65) 


144.291 


13.807 


55° 


1 


7.0 (0.5, 12.0) 


0.55 (0.40, 0.65) 


0.25 (0.20, 0.50) 


133.730 


16.099 




2 


10.0 (8.0, 11.0) 


0.35 (0.35, 0.40) 


0.50 (0.30, 0.50) 


225.629 


37.159 




3 


9.0 (4.0, 11.0) 


0.35 (0.25, 0.40) 


0.45 (0.40, 0.65) 


183.933 


37.418 




4 


5.0 (0.0, 15.0) 


0.65 (0.30, 0.70) 


0.20 (0.20, 0.65) 


144.406 


11.170 



I Sakamoto et al. 1999 1 the ratio would be approximately 1 . This 
is consistent with the ratio of 4 x 10"' obtained using the Mueller 



[Sanchez et aTlpOOg) conversion, considering that both estimates 
have large uncertainties. A more precise way to constrain the cold 
H2 masses is by more direct measurements of CO { J = 1-0) emis- 
sion in the millimeter range. [Heifer et al.| ( [2003^ measured a peak 
molecular surface density (i.e. the peak from the 6 arcsec beam 
size) of 815 Mq pc"^. Over a 3 arcsec aperture a total gas mass 



of approximately 1.5 x 10' Mq would then be expected. Sakamoto 
|erar](fT999| (~ 3 arcsec beam size) report a molecular mass of 
4 X 10** Mq (2.67 X 10^ Mq when using the CO-to-H2 conversion 
factor of [Heifer et al.|2003) within a 15 arcsec diameter aperture. 
Fig.[24|shows the measured gas mass distribution compared to (1) 
the dynamical mass distribution and (2) the gas masses obtained 
from CO (y = 1-0) measurements. We used figure 2 and equa- 
tion 1 of [Sakamoto et al.|jl999[ > with the CO-to-H2 conversion fac- 
tor of [Heifer et al. ( 2003 ^ to estimate the gas mass distribution at 



smaller radii. Within a 3 arcsec diameter aperture the gas masses 
from [Sakamoto et al.[ ljT999| and [Heifer et al.[ ( |2003| l agree very 
well. This shows that we overestimated the cold gas mass traced by 
hot H2 by a factor of ~ 11 when using the conversion of [Muellerj 
[Sanchez et aL] ( [2006} . When comparing the gas mass profile (cal- 
ibrated to match the gas mass at r = 1.5 arcsec derived from CO) 
with the dynamical mass, we find that within a 3 arcsec aperture the 
gas is approximately 5% of the dynamical mass. At smaller radii 
this fraction is even lower, as the gas mass profile is steeper than 
the dynamical mass profile. The gas mass distribution in quadrants 
1 and 4 is larger by a factor ~ 1.5 than in quadrants 2 and 3. How- 
ever, due to the small total mass we do not expect that the irregular 
gas distribution has a large influence on the stellar kinematics and it 
is unlikely the cause for the asymmetry in the stellar kinematics; it 
might, however, be related to whatever is responsible for the latter 

Concerning the results of the black hole mass measurement, 
the gas mass has no significant influence. As the gas mass distribu- 
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Figure 20. as a function of M. (left column, marginalised over Tbuigc and Tjisc), Tbuige (middle column, marginalised over M, and Tdi,sc) and Tjisc (right 
column, marginalised over M, and Tbuige) for NGC 3368. Different inclinations are indicated by different colours (red; i = 52°, black: i = 53°, green: i = 54°, 
blue: i = 55°). 



tion has a similar radial profile as the stellar mass, our T^uigc likely 
includes the gas mass, such that the true stellar mass-to-light ratio 
is maybe slightly lower than our nominal Tbuige- But this is not a 
serious problem for the modelling. However, the fact that the gas is 
not evenly distributed and the clouds already have a mass of order 
10' Mq each, can weaken the evidence for the presence of a black 
hole, as it would imply that the centre is slightly out of equilibrium. 



5 DYNAMICAL MODELLING OF NGC 3489 

5.1 Construction of the Luminosity Profile for Modelling 

Given the apparent similarity of NGC 3489's inner structure to 
that of NGC 3368 (modulo the presence of a secondary bar in 
NGC 3368), including the strong isophotal twist created by the bar 
in NGC 3489, we followed a similar strategy for constructing the 
luminosity profiles. That is, we divide the galaxy into separate disc 
(which includes the discy pseudobulge) and central classical bulge 
components, with the disc treated as having a constant observed el- 
lipticity of 0.41. The disc surface brightness profile is an azimuthal 
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Figure 21. The major axis kinematics of NGC 3368 is sliown as black points for quadrants 1 to 4. Overplotted is the fit of the best model with black hole. 



average with fixed ellipticity down to r = 4.9 arcsec, with the pro- 
file at smaller radii being the extrapolated inner-exponential fit from 
Fig.|9] 

The classical bulge profile is the result of a free-ellipse fit 
to the inner-disc-subtracted WFPC2 F814W image. The latter was 
created by generating a model disc with ellipticity = 0.41 and pro- 
file matching the exponential part of the fit in Fig. [9] (scale length 
= 4.9 arcsec), and then subtracting it from the dust-corrected PC 
image. 

The deprojection was done in the same way as for NGC 3368 
for bulge and disc component separately. 

5.2 Dynamical Models 

NGC 3489 has only a weak large-scale bar and no nuclear bar. The 
measured kinematics and line indices are largely symmetric apart 
from the asymmetry in v. Thus non-axisymmetries are not expected 
to play a role as big as in NGC 3368. We first use only SINFONI 
data to model all four quadrants separately. However, we expect 
that, as for NGC 3368, due to the small field of view of the SIN- 
FONI data it will be difficult to constrain T^isc, as the data cover 
only that part of the galaxy where the classical bulge dominates. 
Thus we try to constrain Tjisc beforehand by modelling SAURON 
and OASIS data alone. As the SAURON data have a large field 
of view including the bar, we use just the inner 10 arcsec for that 
purpose. Finally, we model the combined SINFONI plus OASIS 
and/or SAURON dataset. 

We do not calculate models for different inclinations, as the 
inclination is well determined from the photometry (/ = 55°). As 



shown for NGC 3368 in the previous section, the inclination cannot 
be constrained better via dynamical modelling and the differences 
within a small inclination range of a few degrees are small (see Tab. 
|4]andFig.|20|. 



5.3 The stellar mass-to-light ratio of the disc 

In order to constrain T^isc we first calculate models using sym- 
metrised SAURON and OASIS kinematics separately. We only cal- 
culate inodels with M. = for the SAURON data, but vary M. be- 
tween and 1.3x 10' Mq for the OASIS data. Fig.|25^ shows Kx^ as 
a function of Tbuige and tusc for the SAURON models. Tjisc is well 
constrained, but in T^uigc a very large range between ~ and ~ 0.68 
is possible. This is due to the fact that the classical bulge is only just 
resolved with SAURON {R™ = 1.3 arcsec, SAURON spatial res- 
olution = 1.1 arcsec). The best-fitting model has Tbuije = 0.28 and 
Ydisc = 0.44. For the OASIS models (Fig. |25]j) the resulting Tdisc 
is higher, which could be a result of the higher cr of the OASIS 
data compared to SAURON. Due to the higher spatial resolution 
(0.69 arcsec), Tbuige is better constrained. The best-fitting model 
has Tdisc = 0.6 and Tbuige = 0.36. It is not possible to constrain M, 
with the OASIS data alone (see Fig.|28l. 



5.4 The black hole mass 

To derive the mass of the SMBH we first use the SINFONI kine- 
matics alone to model the four quadrants separately. We chose a few 
values for Tjisc around 0.44. The results with the corresponding 3cr 
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Figure 22. difference between the best-fitting model without black hole and the best-fitting model with black hole (A^^ = 2, = 2?iiCt?noBH ~ 
over all 21 velocity bins) for all LOSVDs of the four quadrants of NGC 3368. Bins where the model with black hole fits the LOSVD better are plotted in green, 
the others in orange. For four bins along the major axis the LOSVDs (open circles with eiTor bars, normalized as in Gebhardt et al.|2000a) and both fits (with 
black hole, full green line and without black hole, dashed orange line) are shown with the corresponding plotted below the LOSVDs. 
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Figure 24. Enclosed mass as a function of radius in NGC 3368 for molec- 
ular mass (black line) and stellar dynamical mass. For the latter, we plot 
values estimated from our modelling of the individual quadrants (Q 1 : red, 
Q2: green, Q3: blue, Q4: cyan). The molecular gas mass is estimated from 
the hot H2 emission in the SINFONI data, converted into a cold gas mass 
using |Mueller Sanchez et ar] j2006 1. The cold gas mass distribution directly 
measured from CO emission by Sakamoto et al. (1999 1 is marked by crosses 
and the cold gas mass from BIMA SONG ^Helfer et al.|2003^ is marked as 
an 'x'. The cold gas mass derived using the hot H2 conversion is clearly 
overestimated compared to the mass derived from CO. The vertical dashed 
line marks approximately the radius of the sphere of influence, defined as 
rsoi = GM./0-2. 



errors are given in Table |5] The mean black hole mass for the four 
quadrants is {M.) = 4.25 x 10* Mq (rms(M.)= 2.05 x 10* Mq). M. 
clearly anticorrelates with Ybuige, but as in the case of NGC 3368 
it does not depend on the specific choice of Tjis^. The mean black 
hole mass for any fixed T^isc is consistent with the result for any 
other Ydisc within Icr. 

The error bars are large, such that a wide range of black hole 
masses is allowed. A solution without black hole is allowed in three 
quadrants within 2 - Act and in one quadrant even within Icr. Thus 
there is no evidence for the presence of a SMBH in one quadrant, 
and only weak evidence in the others, when modelling the SIN- 
FONI data alone. The fit of the best model in each quadrant to v, cr, 
Ht, and /^4 along the major axis is shown in Fig.|26| 

The resulting M,, Tbuigc and Tjisc of the four quadrants agree 
with each other within < 2 cr, and as there are also no obvi- 
ous strong inconsistencies between the kinematics of the quad- 
rants, we fold the LOSVDs of the four quadrants (the LOSVDs 
of quadrants 1 and 4 were also flipped, such that v and change 
sign). For the folded data we find a best-fitting black hole mass of 
M, = 5.0x 10* Mq at Tbuigc = 0.56. This is in good agreement with 
the results of the individual quadrants. A solution without black 
hole is allowed within 3cr, thus as a conservative result we can only 
give a 3cr upper limit of 1.3 X 10^ Mq for the SMBH in NGC 3489, 
when using just the SINFONI data. 

The non-dependence of M. on Tdisc can be explained by the 
very small field of view of the SINFONI data, which covers only 
the very central part of the galaxy, dominated by the classical bulge. 
This might also explain the relatively weak detection of a SMBH 
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Figure 25. ^xl = -X^^\^ as a function of Tjisc and Tbuigc for (a) the sym- 
metrised SAURON data and (b) the symmetrised OASIS data of NGC 3489. 
The coloured regions are the 1 - 5 cr confidence intervals for two degrees 
of freedom. Each calculated model is marked as a black dot. 



in NGC 3489 despite the high quality data. It therefore seems rea- 
sonable to include kinematics at larger radii, like the SAURON 
or the OASIS kinematics, as these datasets cover a large fraction 
of the disc and therefore are able to constrain Ydisc very well, as 
shown above. We should keep in mind however, that the SAURON 
and OASIS velocity dispersions do not fully agree with each other, 
are larger than the SINFONI dispersion and show some deviations 
from axisymmetry, which might possibly be due to the strong dust 
features. In order to determine how strong these differences affect 
the result of the modelling we do three sets of models: the first 
one with SINFONI and OASIS data (using OASIS data between 
r = 0.5 arcsec and 4 arcsec), the second one with SINFONI and 
SAURON data (using SAURON data between r = 1 arcsec and 
10 arcsec) and the third one with all three datasets (with OASIS 
data between r = 0.5 arcsec and 4 arcsec and SAURON data be- 
tween 4 arcsec and 10 arcsec). 

Fig. [27] shows the resulting Ax^ contours for the combined 
SINFONI, OASIS and SAURON data. The error contours are very 
narrow and both Tbuigc and Tjisc are very well constrained. Fig. 
28 shows the resulting Ax^ profiles for all data combinations we 
used. It is clear that the mass of the black hole can be much bet- 
ter constrained when including SAURON and/or OASIS data. The 
constraints on Tbuigc and Tdisc are also much stronger in these cases. 
Using SAURON data in addition to SINFONI and OASIS does not 
seem to improve the measurement of M. and Tbuigc- The scatter 
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Figure 26. The major axis kinematics of NGC 3489 is shown as black points for quadrants 1 to 4, the folded data and (in the last column) the folded SINFONI, 
OASIS and SAURON data. Overplotted is the fit of the best models with black hole. 



in the /S.x^ profiles is quite large for the models of the combined 
datasets, despite the good quality and high S/N of the individual 
datasets and despite the comparatively small scatter in the mod- 
els of individual datasets. The uncertainties of the SMBH mass 
measurement therefore do not seem to be dominated by statistical 
errors, but instead by systematics. Systematic errors can be intro- 
duced e.g. due to the differences in the kinematics of the individual 
datasets. Systematic errors in the modelling (e.g. slightly different 
results for different quadrants) could add to the scatter as well, but 
are difficult to quantify. We measure the formal 1 a errors (corre- 
sponding to l^x^ = 1 for one degree of freedom) by fitting a third 
order polynomial to each curve in Fig. [28] The best values for M., 
Tbuigc and Tjis^ given in Table |6] refer to the minimum of the fit 
and the associated lix^ ^ 1 region. We then determine the sys- 
tematic error, introduced by the differences between the datasets, 
from the scatter of the best fits for models with combined datasets. 
Thus when using all available data, we obtain a black hole mass 
of M. = {(>.W_l\%a ± 0.64|sys) X 10*^ M0, a bulge mass-to-light 



ratio Yi 



0.45 ± 0.02|s,a, ± 0.03|sys and a disc mass-to-light ratio 



I disc ■ 



bulge 



5.5 Evidence for a black hole in NGC 3489 

In order to illustrate where the influence of the black hole is largest. 
Fig. 29 shows the^- difl"erence between the best-fitting model with- 
out a black hole and the best-fitting model with a black hole for all 
LOSVDs of the combined SINFONI+OASIS+SAURON dataset. 
The fit with black hole is generally better in 80% of all bins. Along 
the major axis the largest differences appear in the LOSVD 
wings, both at negative and positive velocities. Improvements of 
the fit appear at all radii. 

Fig. 30 shows hx^ summed over all angles and velocities as 



a function of radius. The l^x^ increase is steepest in the region 
covered by the SINFONI data (A;('5jfjpQj,jj ~ 72, mostly coming 
from the region < 0.15 arcsec, corresponding to ~ 3rsoi). It then 
grows by about the same amount in the region of the OASIS data 
(0.5 - 4 arcsec). At larger radii (region of the SAURON data) the 
increase is only small. 



Although with the OASIS data alone it is not possible to con- 
strain the mass of the SMBH, the region covered by this dataset 
(~ 0.5 - 4 arcsec) seems to be crucial for the lower limit on M., 
which is not possible to derive with the SINFONI data alone. This 
means that differences between models without black hole and 
models with black hole (say, M, = 6 x 10* Mq) should not only 
appear within the sphere of influence, but also further outside. This 
should not be surprising. For example, some of the effects of a cen- 
tral mass concentration in an isotropic system can be mimicked in a 
system without such a concentration by enhanced radial anisotropy 
(e.g., |Binney & Mamon|1982} . In such a case, the region where ra- 
dial anistropy is required extends well outside the nominal sphere 
of influence of the central mass. Only in cases where the best-fit 
model without a black hole has exactly the same orbital structure 
as the best-fit model with the black hole would differences between 
the fits be (roughly) confined to the sphere of influence. This is 
in agreement with the observations in NGC 3368, Fornax A and 
NGC 4486a jNowak et ar|2 007 , 2008 1, where a general improve- 
ment of the fit even outside the sphere of influence, was observed. 



The total A;^'^, summed over all LOSVDs, between the best- 
fitting model without a black hole and the best-fitting model with a 
black hole is given in the last column of Table [6] 
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Figure 28. A;(^q = x ~ xf^^i^ ^ function of M, (left column, marginalised over Ttuije and Tjisc), Tbuige (middle column, marginalised over M. and T^n^c) 
and Tdisc (right column, marginalised over M, and Tbuiee) for NGC 3489. The colours indicate the datasets used for the modelling (black: SINFONI, grey; 
OASIS, blue: SINFONI+SAURON, green: SINFONI+OASIS, red: SINFONI+OASIS+SAURON). 
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Figure 29. Same as Fig.|22]for NGC 3489 for the models using the averaged SINFONI+OASIS+SAURON data. Instead of a linear radial scale we used jc''^ 
and due to the increasing bin size with radius. The conversion is shown in the bottom scale. 



5.6 Discussion 

No attempts have been made in the literature to measure the mass of 
the SM BH in NGC 3489. From the M.-ar relation of |Tremaine et al.| 



1 2002 I we would expect a mass between M. 



B.9 km s derived from the SAURON data ffimsellemetal. 



2004 1 and M. = 9.2 x lO*" Mq for cr^ = 102.5 km s ' derived from 



the OASIS data jMcDermi d et al.'2006^. From the relation between 
M. and A'-band magnitude ( Marconi & Hunt,2003^ we would ex- 



pect a black-hole mass of M. = 1.92 x 10' Mg if it correlates 
with the total (photometric) bulge magnitude ^K,ioVd\ ~ 

-21.91, 

or 4.94 X 10* Mq if it correlates with the classical bulge magnitude 
Mjf.buige = -20.60 only. 

The stellar population models of 'Mar aston| ( [1998] [200 5^ sug- 
gest an age of ~ 1 Gyr for the best-fitting Tbuige = 0.45 ~ Ydisc and 
a high metallicity [Z/H]~ 0.67, and an age of ~ 2-3 Gyr for a solar 
metallicity population (assuming a Salpeter IMF). This is in agree- 
ment with |McDermid et al.| ( |2006) , who find both an age gradient 
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Figure 27. Same as Fig. [19] for NGC3489, with an inclination i = 55°. 
Tlie averaged SINFONI data, OASIS data between 0.5 and 4 arcsec and 
SAURON data between 4 and 10 arcsec were used for the modelling. 
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Figure 30. Same as Fig.|23]for the combined SINFONI+OASIS+SAURON 
data of NGC 3489. 



(from ~ 2 - 3 Gyr in the outer regions to ~ 1 Gyr in the centre) and 
a metallicity gradient (from ~solar in the outer regions to ~ 0.6 in 
the centre). It is also compatible with |Sarzi et al.|j2005 ), who found 
a mean age of ~ 3 Gyr assuming solar metallicity. The stellar mass 
within rsoi ~ 0.053 arcsec is a 6 x 10' Mg. If the best-fitting mass 
for the black hole of M, = 6.0 x 10* Mg were entirely composed of 
stars, the mass-to-light ratio would increase to 0.9. This would be 



typical for an older stellar population (~ 6 Gyr for a high metallic- 
ity [Z/H]=0.67 and a Salpeter IMF), and therefore conflict with the 
values found by |Sarzi et aLl ( |2005) and |McDermid et al.| ( (2006l l. 



6 SUMMARY AND DISCUSSION 

We analysed near-IR integral-field data for two barred galaxies that 
host both a pseudobulge and a classical bulge component. Both 
galaxies show fast and regular rotation and a cr-drop at the cen- 
tre, which in the case of NGC 3368 is more pronounced and may 
have developed from gas, transported to the inner region by the 
bars and spiral arms. The kinematics of NGC 3368 - in particular 
the velocity dispersion - is asymmetric. The reasons for that could 
be dust or (less likely) the non-axisymmetric potential induced by 
the two bars. The gas distribution is also inhomogeneous, but as 
the total gas mass accounts for only < 5% of the dynamical mass, 
this has probably no significant influence on the stellar kinematics. 
There are two kinematically decoupled gas clouds located a few 
tens of parsecs north of the centre. Each cloud has a total mass of 
order 10' Mq. The stellar kinematics of NGC 3489 is very regular, 
with a slight asymmetry in the velocity field. All other kinematic 
parameters and the line indices are consistent with axisymmetry. 
No gas emission was found in NGC 3489. The near-IR line indices 
Nal and CO show a negative gradient in both galaxies, indicating 
an age and/or metallicity gradient. 

We applied axisymmetric dynamical models to derive the 
SMBH masses in NGC 3368 and NGC 3489. In our models we as- 
sume that the galaxy potential can be decomposed into three com- 
ponents: the central black hole, an inner, classical bulge (with mass- 
to-light ratio Tjjuige) and the disc (T^isc); the disc component in- 
cludes the pseudobulge. The inclination of the models is fixed by 
the isophotes of the outer disc. For NGC 3368 we modelled the four 
quadrants of our IFU data independently and the resulting black 
hole masses and mass-to-light ratios agree very well. We find that 
M. is largely independent of Ynsi: and anticorrelates with T(,uigc- 
The average black hole mass for the four quadrants and an inclina- 
tion (■ = 53° is (M.) = 7.5 X 10' Mq (rms(M.)= 1.5 x 10' Mq). A 
solution without a black hole is excluded by a; 4-5(7". The errors, 
however, cover a large range in M.. The largest uncertainty for M. 
comes from the unknown Tbuigo and independent constraints, e.g. 
from stellar population modelling, would likely improve the results. 
However, unless the shape of the IMF is known, mass-to-light ra- 
tios from stellar population analyses are ambiguous. The scatter 
from quadrant to quadrant is smaller than the uncertainty related to 
Tbuigc, suggesting that the symmetry assumption plays a minor role 
for the uncertainty of M.. Our results do not significantly depend 
on the inclination (within the photometrically allowed inclination 
range). 

For NGC 3489, modelling of the four SINFONI quadrants 
likewise gave consistent black hole masses and mass-to-light ra- 
tios. Similar to NGC 3368 the errors in M, are large, the black 
hole mass is independent of Tjisc and it clearly anticorrelates with 
Tbuigc. Modelling the folded SINFONI data gives the same result 
as for the individual quadrants; thus, non-axisymmetries do not 
seem to play a role. When including OASIS and/or SAURON data, 
Tbuigc and therefore also M, could be much better constrained. Us- 
ing all three datasets, we derived for NGC 3489 a SMBH mass of 
M. = {6Wt_oT4Ua ± 0.64|sys) X 10' Mo with a bulge mass-to-hght 
ratio of Tbuigc = 0.45 ± 0.02|s,a, ± 0.03|sys and a disc mass-to-light 
ratio Td,sc = 0.47+55[5^|s,a, ± 0.05|sys. A solution without a black hole 
is excluded with high significance. To derive a firm lower limit to 
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Table 5. Resulting black hole masses M, and //-band mass-to-light ratios T|,ui„j. and Tjis^ of NGC 3489. The lower and upper 
3o" limits are given in brackets. The total of the best model with black hole and the;j'^ difference between the best model 
without black hole and the best model with black hole are given in the last two columns. 



Quadrant 


M. [10'^ Mq] 


bulge 


Tdisc 


2 

^ mill 




1 


4.0 (0.0, 8.0) 


0.60 (0.44, 0.76) 


0.36 (0.28, 0.64) 


106.812 


3.289 


2 


1.0(0.0,7.0) 


0.60 (0.44, 0.80) 


0.60 (0.28, 0.64) 


100.554 


0.301 


3 


6.0 (0.0, 13.0) 


0.48 (0.28, 0.72) 


0.52 (0.28, 0.64) 


57.612 


10.487 


4 


6.0(1.0, 10.0) 


0.52 (0.36, 0.64) 


0.32 (0.28, 0.64) 


80.077 


16.786 


folded 


5.0 (0.0, 13.0) 


0.56 (0.28, 0.72) 


0.52 (0.28, 0.64) 


47.877 


9.060 



Table 6. Resulting black hole masses M, and //-band mass-to-light ratios Tbuigc and Tjisc of NGC 3489 for the folded 
SINFONI data alone and in combination with SAURON and OASIS kinematics. The low er and upper lo" limits (1 degree of 
freedom), determined by fitting a third order polynomial to the profiles of Fig. |28| are given in brackets. The total 
of the best model with black hole and the difference between the best model without black hole and the best model with 
black hole are given in the last two columns. 







M. [10* Mq] 






2 




SINFONI 




5.97 (3.64, 8.13) 


0.53 (0.48, 0.58) 


0.55 (0.48, 0.62) 


106.812 


3.289 


SINFONI + SAURON(l 


- 10") 


4.56(4.03,5.12) 


0.52 (0.50, 0.54) 


0.36 (0.32, 0.40) 


368.458 


81.83 


SINFONI + OASIS(0.5 - 


4") 


5.81 (5.21,6.46) 


0.46 (0.44, 0.48) 


0.52 (0.50, 0.55) 


551.297 


139.86 


SINFONI + OASIS(0.5 - 


4") + SAURON (4 - 10") 


6.00 (5.46, 6.56) 


0.45 (0.43, 0.47) 


0.47 (0.45, 0.48) 


606.653 


173.86 



M., data between ~ 0.5-4 arcsec seem to be crucial, in addition to 
the high-resolution SINFONI data in the centre. With OASIS data 
alone, no limits on M. could be placed. There are some inconsis- 
tencies in the kinematics between the three datasets, which seem to 
be the main source of systematic errors. In particular when mod- 
elling OASIS data alone, we get a higher Tjisc than if modelling 
SAURON data alone (because the inner cr is higher in the OASIS 
data than in the SAURON data). 

The implications for the M.-cr relation and the M.-M^ rela- 
tion are illustrated in Fig.[3T] For NGC 3368 the mean M. of the 
four quadrants and the rms, and for NGC 3489 M, from the com- 
bination of SINFONI, SAURON and OASIS data with its statis- 
tical Icr error is plotted against cr and Mk using the relations of 
[Tremaine et~arHjo02l >, [Ferrarese & Fordl ( [2005| >, |Marconi & Hunt] 



and 



Graham 



( |2007 l|. 



All values for cr. and ct./r were mea- 



([2003) 

sured using the effective radius of the total photometric bulge, as 
was done for all the galaxies contributing to the [Tremaine et aL] 
( |2002^ and [Ferrarese & Ford| ( |200?| relations. No attempt to de- 
termine (Te for the classical components has therefore been made, 
but as we use luminosity-weighted measurements, all values de- 
termined from high-resolution data represent mostly the classical 
bulge. 

The agreement of NGC 3368 with the M.-cr relation largely 
depends on the value of cr which is used. The small cr = 
98.5 km s"' measured within the SINFONI field of view is in 
good agreement with the M.-cr relation. When combining the SIN- 
FONI cr with cr measurements of Whitmore, Schechter & Kirsh-| 
[ner| ( (T979) , rHeraudeau et al.|(T999| and |Vega Beltran et al.|pOOT) 
a v£ilue of cr^ii = 104 km s" is obtained. The velocity disper- 
sions from t he literature alone however (e.g. cr^ = 117 km s ' 
estimated by Sarzi et ar]|2002 cr^ = 130.9 km s~' respectiv ely 
cr,/8 = 129.9 kms ' 
cr s= 150 km s"' by 



measured b' 



Moiseev 



red by Her; 
et al.|200¥ 



Heraudeau et al. 



1999 



are significantly larger 
than expected by this estimate and not or only marginally in agree- 
ment with the M.-cr relation. With a K-hanA magnitude of -23.42 
for the total photometric bulge, NGC 3368 falls far (a factor of 
~ 12) below the M.-Mk relation of [Marconi & Hunt] (2003) . If 



we postulate that the SMBH only correlates with the magnitude of 
the classical bulge, the situation improves. With M™ = -19.48 
NGC 3368 now lies a factor of ~ 5 above the M.-M^ relation of 
[Marconi & Hunt 1 2003|l, b ut is in good agreement with the M.-Mk 
relation of [Graham[p007l >. 

For NGC 3489 the situation is similar. M. is in excellent agree- 
ment with the M.-cr relation when using either the SINFONI mean 
cr = 91. 1 km s"' or the SAURON values cr^ = 88.9 km s"' and 
cre/8 = 94 km s"' . It is still in reasonably good agreement with the 
relation when using the OASIS measurements (cr^ = 102.5 km s"', 
CTc/g = 108.9 km s"') or when taking into account other cr measure- 
ments from the literature (cr^/g = 115 km s"' using Whitmore et al. 



[I979|[Dalle Ore et al.[199Tl[Smith et al.[2000[[Barth et al.[2002[ and 
the SINFONI value). With a K-hmd magnitude of the total photo- 
metric bulge of M™ = -21.91 NGC 3489 also falls far below the 
M.-Mk relation of [Marconi & Hunt| ( |2003) and [Graham] pOOT) !, 
but is in excellent agreement if the magnitude of the classical bulge 
component is considered (Mk = -20.60). 

The large difference in the cr measurements makes it diffi- 
cult to draw any firm conclusion with respect to the location of 
pseudobulges in the M.-cr relation, and at the same time illustrates 
that measurement errors in cr may play a larger role than one may 
have thought, in particular when dealing with small galaxy sam- 
ples. NGC 3368 would fall far below the M.-cr relation when opti- 
cal longslit kinematics alone are used. These discrepancies between 
the cr measurements might at least partly be due to dust, which af- 
fects the optical data much more than the near-IR data. 

The A'-band magnitudes on the other hand can be determined 
very accurately even for subcomponents of the galaxy. Taken at 
face value, both galaxies clearly do not follow the M.-Mk relation 
of [Marconi & Hunt|j2003^ when considering the /f-band magni- 
tudes of the total photometric bulge, but are in better (NGC 3368) 
or even excellent (NGC 3489) agreement with it when considering 
the classical bulge magnitude only. 

If we take into account that a stellar population becomes 
fainter when it ages passively (2.3 mag in K band for a solar metal- 
licity population between 1 and 10 Gyr, based on the stellar pop- 
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ulation models of |Maraston|1998||2005^ , the pseudobulges would 
move toward the M.-Mk relation with time. Given the uncertainties 
on the age estimate, the exact size of the effect is unclear. Keeping 
in mind this caveat, this is in line with Greene et al.| ( |200 8X who 
conclude that pseudobulges follow the M.-cr relation, but not the 
M.-M(,iiige relation, as well as with |Gadotti & Kauffmann] ( |2009[ l, 
who find that pseudobulges follow only one of the two relations, if 
any. In order to strengthen our results, studies of a larger sample of 
pseudobulges similar in design are necessary. 

Whether modelling single quadrants of obviously non- 
axisymmetric galaxies with an axisymmetric code is a good ap- 
proximation and gives the correct black hole masses is certainly 
still an issue that remains to be resolved. The recently developed 
triaxial codes of |de Lorenzi et ar|p007[ > and |van den Bosch et al.| 
( |2008| l will have the potential to solve this issue in the future. 
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Figure 31. Left panel: Location of NGC 3368 and NGC 3489 with respect to the M.-cr relation (black: |Tremaine et al.|2002| red: [Ferrar ese & Ford|2005^. 
The velocity dispersion from our SINFONI measurements (open symbols) and values for ctj (filled black symbols) and CTe/g (filled red symbols) derived from 
the literature are plotted for each galaxy. Right panel: Location of the two galaxies with respect to the M.-Mg relation of [Marconi & Hunt^ (2003 1 (black) 
and [GrahamH200'7) (red). The ^T-band magnitudes of the total photometric bulges are plotted as filled black symbols, the magnitudes of the classical bulge 
components as filled red symbols. In both panels the average M. of the four quadrants and the rms is plotted for NGC 3368, and M. from the combination of 
SINFONL SAURON and OASIS data with its statistical Icr error is plotted for NGC 3489. 
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